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ABSTRACT 


The  gas  dynamic  structure  of  high  Edtitude  rocket  plumes  has  been 
studied  in  a  low  density  tunnel  to  obtain  basic  data  on  plume  size  and 
shape  and  to  determine  the  extent  of  the  mixing  region  between  the  free- 
stream  and  rocket  gases.  The  rocket  was  simulated  by  specially  de¬ 
signed  ejectors  using  helium  and  CO2  as  the  gases,  with  conditions 
chosen  to  simulate  a  100,  000 -lb-thrust  rocket  at  an  altitude  of  80  km. 
Data  taken  included  impact  pressure  surveys,  density  measurements 
of  the  free-stream  and  plume  gas,  temperature  measurements,  and 
color  photographs.  The  last  three  of  these  were  taken  using  electron 
beam  techniques.  The  result  of  the  study  is  a  set  of  basic  data. 
Analysis  and  interpretation  of  the  data  will  be  performed  under  a  sepa¬ 
rate  study. 


AeDC-TR-71.2S 


CONTENTS 

Page 

ABSTRACT .  iii 

NOMENCLATURE . viii 

I.  INTRODUCTION .  1 

U.  BASIC  THEORY  AND  SELECTION  OF  TEST  CONDITIONS 

2. 1  General  Featiires  of  a  Htmie  .  . .  2 

2.  2  Simplified  Descriptions  of  the  Flow  Field .  4 

2.3  Selection  of  Test  Conditions .  14 

in.  DESCRIPTION  OP  EQUIPMENT 

3.  1  Tunnel  M .  17 

3.2  T unnel  D .  25 

3.3  Plume  Generation  System .  25 

3.4  Tunnel  Conditions  and  Pressure  Measuring  Equip¬ 
ment  .  31 

3.  5  Electron  Beam  Apparatus .  33 

3. 6  Photographic  Apparatus .  413 

IV,  DATA  ACQUISITION  AND  REDUCTION  TECHNIQUES 

4, 1  Tunnel  Conditions- .  46 

4,  2  Plume  Sowce  Conditions .  49 

4.  3  PLume  Impact  Pressure  Surveys .  49 

4. 4  Electron  Beam  . .  51 

V.  EXPERIMENTAL  RESULTS 

5.  1  Photographs .  55 

5.  2  Impact  Pressure  Surveys .  60 

5.3  Density  and  Temperature  Measurements .  90 

5.4  Comparison  of  Results . 115 

VI.  SUMMARY  AND  CONCLUSIONS .  125 

REFERENCES .  126 

ILLUSTRATIONS 

Figure 

1.  General  Features  of  a  High  Altitude  Rocket  Plume ...  3 

2.  Flow  Field  for  Expansion  into  a  Region  of  Uniform 

Pressure .  5 

3.  Plume  Spreading  Parameter,  A .  8 

4.  Simplified  Sketch  of  Plirnie  Contact  Surface .  9 


V 


AEDC-TR-71-25 


Figure  Page 

5.  Approximation  Representation  for  1  “COS  0 .  13 

6.  Characteristics  Solutions .  19 

7.  Tunnel  M  .  . .  23 

8.  Tunnel  D .  26 

9.  Photograph  of  Plume  Generators  .  27 

10.  Sim\d.ator  and  Impact  Pressure  Probe  Installation  ...  28 

11.  Pliune  Gas  Supply  System .  29 

12.  Instrumentation  Schematic .  32 

13.  Tunnel  D  Simulator  and  Impact  Pressure  Probes  during 

Typical  Run .  34 

14.  Photograph  of  Electron  Beam  Apparatus  in  Tunnel  M.  .  36 

15.  Electron  Gun  Power  Supply  Schematic .  37 

16.  Electron  Beam  Schematics .  38 

17.  Photograph  of  Current  Collector .  41 

18.  Photograph  of  Spectrometer .  42 

19.  Downstream  View  of  Electron  Beam  -  Spectrometer- 

Telescope  Setup .  43 

20.  Integrator  Schematic .  45 

21.  Transmission  Characteristics  of  Filters .  47 

22.  Typical  Rotational  Spectra .  52 

23.  Boltzmann  Hot,  DWM  0020  .  54 

24.  Photograph  of  Hume  for  Condition  ICS,  Tunnel  M, 

He  Hume,  Aq/A*^*  =  1.  61,  d""  =  0.  315  in. ,  and 

Pc  =  5  psia .  57 

25.  Photograph  of  Hume  for  Condition  IC2,  Tunnel  M, 

He  Hume,  Ag/A'‘‘  =  1..61,  d'*'  =  0.  315  in. ,  and 

p^  =  2  psia .  57 

26.  Photograph  of  Hume  for  Condition  IB,  Tunnel  M, 

-He  Hume,  .A^A'*'  f_25,  d"^  =  0.030  in. ,  and 

Pc  =  450  psia .  58  “ 

27.  Photograph  of  Hume  for  Condition  lAB,  Tunnel  M, 

COg  Hume,  A^/A*  =  25,  d*  =  0.  030  in. ,  and 

Pc  =  400  psia . ■ .  58 


VI 


AEDC-TR-71-25 


Figure  Page 

28.  Photograph  of  Argon  Plume,  Tunnel  M,  A^/A*  =  1.  6, 

d*  =  0.  315,  and  =  3  psia .  59 

29.  Impact  Pressure  Starveys  for  Condition  IC5 .  61 

30.  Impact  Pressure  Surveys  for  Condition  lAB .  65 

31.  Impact  Pressure  Surveys  for  Condition  IIB .  70 

32.  Impact  Pressxire  Surveys  for  Condition  IC2 .  73 

33.  Impact  Pressure  Surveys  for  Condition  lAC .  76 

34.  Impact  Pressure  Surveys  for  Condition  lAA .  81 

35.  Impact  Pressure  Surveys  for  Condition  IB .  86 

36.  Electron  Beam  Number  Densities  for  Condition  ICS  .  .  91 

37.  Electron  Beam  Number  Densities  for  Condition  IC2  .  .  94 

38.  Electron  Beam  Number  Densities  for  Condition  IB 

(x  *=  3  in. ) .  97 

39.  Electron  Beam  Number  Densities  with  Collisional 

Quenching  Effects  for  Condition  ICS .  90 

40.  Electron  Beam  Number  Densities  with  Collisional 

Quenching  Effects  for  Condition  IC2 . 102 

41.  Electron  Beam  Number  Densities  with  Collisional 

Quenching  Effects  for  Condition  IB  (x  =  3  in. )  .  .  .  .  .  105 

42.  Locations  of  the  Maxima  of  the  He  and  N2  Number 

Densities . 106 

43.  Centerline  Number  Density  of  He  versus  Inverse  Square 

of  Axial  Distance . 108 

44.  Rotational  Temperature  of  N2  as  a  Function  of  the 

Number  of  Spectral  Lines  Used . 110 

45.  N2  Rotational  Temperatime  versus  Radial  Position 

Using  Ten  Spectral  Lines . 113 

46.  N2  Vibrational  Temperature  versus  Radial  Position 

for  Condition  IC2 . 116 

47.  N2  Vibrational  Temperature  versus  Radial  Position 

for  Condition  ICS . 116 

48.  Comparison  of  Results  for  Condition  ICS  .  .  .  .-.  .  118 

vii 


AEDC-TR-71-25 


Figure  Page 

49.  Comparison  of  Results  for  Condition  IC2 .  122 

50.  Comparison  of  Results  for  Condition  IB  {x  =  3  in. )  .  .  125 

51.  Normalized  Location  of  Maximum  Impact  Pressure.  .  126 

TABLES 

I.  Summary  of  Test  Conditions  . .  15 

n.  Nominal  Tunnel  M  Nozzle  18  N2  Operating  Condi¬ 
tions  . . .  24 

III.  Plume  Simulator  Dimensions  . .  30 

IV.  T est  Summary . .  50 

V.  Experimental  Data  on  Mach  Disk  Location  as  Deter¬ 
mined  fay  Impact  Pressure  Surveys .  89 

VI.  Electron  Beam  Results .  107 

VII.  Ten  Spectral  Line  Tj^  as  a  Function  of  Radial 

Position  . .  114 

Vm.  Vifarational  Temperature  as  a  Function  of  Radial 

Position .  117 

APPENDIXES 

I.  Correction  of  Pressure  Measurements  for  Orifice  and 

Probe  Effects .  ^31 

n.  Flow  Diagnostics  Using  Excitation  by  an  Electron 

Beam . 

NOMENCLATURE 

First 

Appearance  or 
Definition 


A 


Einstein's  spontaneous  emission  , 
coefficient 


Eq.  (n-1) 


AEDC-TR-71-25 

A* 

Throat  area  of  nozzle 

Eq.  (3) 

-A-e 

Exit  area  of  nozzle 

Fig.  3 

Bo 

Eotational  constant  for  the  v'"  =  0 
state  of  N2X'E 

Eq.  (II-5) 

C,  Cq,  Cj,  C2 

Constants 

Appendix  II 

Cp 

Thrust  coefficient  of  rocket 

Eq.  (5) 

^  m 

Maximum  thrust  coefficient  of  a 
rocket  (exhausting  into  a  vacuum) 

Eq.  (6) 

c 

Speed  of  light 

Eq.  (n-1) 

d* 

Diameter  of  nozzle  throat 

Table  I 

Diameter  of  nozzle  exit  =  2yg 

Table  I 

F 

Vacuum  thrust  of  rocket 

Eq.  (4) 

Maximum  vacuum  thrust  of  rocket 

(Ag/A*-.a.) 

Eq.  (6) 

[Glo 

Muntz  G  “factor 

Eq.  (n-6) 

h 

Planck's  constant 

Eq.  (n-1) 

h 

0 

Stagnation  enthalpy  of  tunnel  gas 

Table  U 

I 

Eadiation  intensity 

Eq.  (II-l) 

1° 

00 

Value  of  Iqq  with  no  coUisonal  decay 
losses 

Eq.  (11-14) 

J— 1 

<  < 

1 

h-* 

Eadiation  intensity  of  vibration- 
rotation  emission 

Eq.  (II-5) 

v' 

^y”'  Vv" 

Eadiation  intensity  of  vibrational 
band  v'v" 

Eq.  (II- 7) 

j 

Current  density  of  electron  beam 

Eq.  (n-1) 

Kl 

Constant 

Eq.  (1) 

K',  K",  K"' 

Eotational  quantum  numbers 

Eq.  (II-4) 

Kn 

Knudsen  number  based  upon  orifice 
diameter,  T^.,  and  wall  conditions 

Appendix  I 

k 

Boltzmann's  constant 

Eq.  (n-4) 

k 

Quenching  constant  for  helium 
intensities 

Eq.  (n-3) 

^0^  ^1 

Quenching  constants  for  v  =  0  and  1 
for  vibrational  band  intensities 

Eq.  (n-14) 

ix 


AEDC-TR-71-25 


L 

Inviscid  scaling  parameter  for  pLume 
geometry 

Eq.  (17) 

Me 

Exit  Mach  number  of  rocket 

Eq.  (9) 

Free -stream  Mach  number 

Eq.  (12) 

m 

Mass  flow  rate  through  rocket  nozzle 

Eq.  (3) 

N 

Nimiber  density 

Eq,  (II- 1) 

P 

Pressure 

p(vv') 

Vibrational  band  strength  for  emission 

Eq.  (n-7) 

p(w') 

Vibrational  band  strength  for  excitation 

Eq.  (n-7) 

p 

Orifice  effect  parameter 

Eq.  (1-3) 

Pc 

Chamber  pressure  of  rocket 

Eq.  (5) 

p.  f. 

Partition  function 

Eq.  (11-5) 

Pi 

Pressure  inside  orifice 

Eq.  (1-3) 

Pifm 

Free  molecular  limit  of  p^ 

Eq.  (1-3) 

Po 

Reservoir  pressure 

Table  H, 
Eq.  (15) 

Po 

Measured  impact  (pitot  probe) 
pressure 

Table  n 

Pw 

Wall  pressure  in  the  absence  of  an 
orifice 

Eq.  (I-l) 

Pi 

Pressure  of  rocket  exhaust  gas  in 
plume 

Eq.  (10) 

Pa, 

Free-stream  static  pressure 

Eq.  (10) 

q(ij) 

Franck-Condon  factor 

Eq.  (n-10) 

Qo, 

Free-stream  dynamic  pressure 

Eq.  (15) 

q 

Heat -transfer  rate 

Eq,  (1-7) 

Free  molecular  value  of  q 

Eq,  (1-9) 

^  FR 

Continuum  (Fay-Riddell)  value  of  q 

Eq.  (I-IO) 

R 

Gas  constant 

Eq.  (1-7) 

R 

Density  dependent  factor  for  quenching- 

Eq.  (11-20) 

corrections  on  vibrational  intensity 
ratios 


X 


AEDC-TR-71-25 


Re 

Electron  transition  moment 

Eq.  (II- 9) 

Reo 

Unit  Reynolds  number  at  a  stagnation 
point  P„V„/a^o 

Table  11 

R®2 

Unit  Reynolds  number  behind  a  standing 
normal  shock  Pa,U„/P2 

Table  II 

rp 

Probe  Reynolds  number  Pa>^a^-pll^2 

Eq.  (I-l) 

Unit  Rejmolds  number  based  on  free- 
stream  condition  PoU^/Ai^j 

Table  I 

r 

Radius  from  or  gin  to  a  point  in  the 
plume  flow  field 

Eq.  (1), 
Fig.  2 

Throat  radius  of  rocket  nozzle 

Eq.  (12) 

r- centroid 

Eq.  (n-12) 

Radius  of  impact  probe 

Eq.  (I-l) 

Tr 

Rotational  temperature 

Eq.  (U-5) 

Tr 

Reference  temperature 

Eq.  (1-6) 

T 

*v 

Vibrational  temperature 

Eq.  (n-7) 

T 

w 

Wall  temperature 

Eq.  (1-2) 

T 

■^0 

Reservoir  temperature 

Table  II 

T 

^op 

Reservoir  temperature  of  plume  gas 

Fig.  12a 

T„ 

Tunnel  free-stream  static  temperature 

Table  11 

Ui 

Velocity  of  rocket  exhaust  gas  in 
plume 

Eq.  (10) 

u„ 

Free-stream  velocity 

Eq.  (10) 

Ulizn 

Limiting  velocity  of  rocket  exhaust 
gas 

Eq.  (4) 

u* 

Velocity  of  gas  at  rocket  throat 

Eq.  (3) 

Up,  ue 

Components  of  velocity  in  plume 

Eq.  (1) 

Fig.  2 

V 

Vibrational  quantum  number 

Appendix  II 

W 

Normalized  distance  for  plume 

Eq.  (12) 

X 

Axial  distance  from  exit  plane  of 
rocket  nozzle 

Fig.  4 

y 

Perpendicular  distance  from  axis  to 
point  in  plume  flow  field 

Fig,  2 

Xi 


AEDC-TR-71-25 


ye 


Value  of  y  for  rocket  ejdiaust 


ymax 

z 

a 


Maximum  size  of  plume 
Parameter  for  orifice  effect 
Thermal  accommodation  coefficient 
Inclination  angle  of  contact  surface 


7  Eatio  of  specific  heats 

6  Boundary-layer  displacement  thickness 

6  Angle  from  axis  to  point  in  the  flow 

field 

Value  of  0  at  edge  of  plume  gas 

Turning  angle  of  rocket  flow  at 
corner 


K 


X 

X 

X 


00 


A 


v'K" 

v"k" 


Nozzle  angle 

Dimensionless  heat-transfer  rate 
parameter 

Parameter  for  nozzle 
W  avelength 

Mean  free  path  at  free-stream  con¬ 
ditions 

Wavelength  of  the  vibration- rotation 
emission  line 


^o 

^2 


V 

p 

Pi 


Coefficient  of  viscosity  at  stagnation 
point 

Coefficient  of  viscosity  behind  stand¬ 
ing  normal  shock 

Coefficient  of  viscosity  at  free-stream 
conditions 

Prandtl-Meyer  angle 

Density  of  rocket  exhaust  at  nozzle 
throat 

Density  of  rocket  exhaust  gas  in  plume 


xii 


Eq.  (14), 

Fig.  4 

Table  I 

Eq.  (1-4) 

Eq.  (1-5) 

Eq.  (10), 

Fig.  4 

Eq.  (8) 

Section  3. 3 

Eq.  (1), 

Fig.  2 

Section  2.  2 

Fig.  4 

Fig.  2 

Eqs.  (1-5), 
(1-7) 

Eqs.  (1),  (7) 

Eq.  (n-1) 

Tables  I 
and  II 

Eq.  (n-5) 

Eq.  (I-IO) 

Table  II 

Table  H 

Eq.  (15) 

Eq.  (3) 

Eq.  (1)  -  - 


AeDC-TR-71-2S 


P2 

Density  behind  a  standing  normal  shock 

Eq. 

(I-l) 

P® 

Free -stream  density 

Eq. 

(10) 

'»'r 

Radiation  life  time 

Eq. 

(n-1) 

A6DC-TR-71-2S 


SECTION  I 
INTRODUCTION 


Thus  far,  attempts  to  account  for  the  observed  radiation  of  rocket 
plumes  at  altitude  have  been  primarily  concerned  with  the  radiative 
characteristics  of  the  ejdiaust  gas  itself,  and  the  role  of  fluid  mechan¬ 
ical  interactions  has  not  been  widely  studied  experimentally.  As  yet, 
a  satisfactory  accounting  for  all  of  the  observed  features  has  not  been 
obtained. 

It  is  possible  at  high  altitudes  that  additional  radiation  can  be  gene¬ 
rated  by  the  interaction  of  the  plume  gas  with  the  free -stream  gas  by 
one  of  several  processes.  The  work  reported  herein  was  undertaken 
to  provide  the  information  necessary  to  arrive  at  a  gas  dynamic  de¬ 
scription  of  the  flow  field  at  high  altitudes.  This  information  can  then 
be  used  with  knowledge  of  the  reaction  rates  and  radiative  character¬ 
istics  for  the  appropriate  components  of  the  gases  to  calculate  the 
radiative  intensity  in  the  interaction  or  mixing  zone  between  the  two 
gases. 

At  lower  altitudes  the  exhaust  plume,  which  is  the  same  order  of 
magnitude  in  size  as  the  rocket,  behaves  as  an  inviscid  solid  body  as 
far  as  the  external  flow  is  concerned.  In  this  regime  the  external  flow 
of  the  atmospheric  gas  and  the  internal  flow  of  rocket  exhaust  gas  may 
be  considered  independently  with  the  boundary  condition  that  at  the 
interface  between  the  two  regions  there  is  equoility  of  pressure  in  the 
two  gas  streams.  The  method  of  characteristics  has  been  used  to 
solve  the  basic  inviscid  flow  equations  for  both  the  external  and  the 
internal  flow  fields.  One  of  the  most  general  of  these  method-of- 
characteristics  programs  was  developed  by  Lockheed  and  is  reported 
in  Ref.  1. 

At  higher  altitudes  the  ambient  flow  becomes  more  rarefied  and 
the  plume  becomes  larger  until  it  far  exceeds  the  rocket  itself  in  size. 
At  this  altitude  the  rocket  plume,  the  shape  and  characteristics  of 
which  are  deteraiined  by  the  rocket  motor  and  properties  of  the  exhaust 
gas  as  well  as  rocket  speed  and  atmospheric  properties,  is  the  major 
factor  in  the  determination  of  the  flow  field  about  the  rocket  and  the 
disturbance  in  the  atmosphere.  "Observables"  such  as  radiation  caused 
by  the  passage  of  the  rocket  through  the  atmosphere  are  largely  deter¬ 
mined  by  the  properties  of  the  rocket  exhaust  plume. 
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At  lower  altitudes  the  radiation  is  primarily  due  to  the  hot  exhaust 
gases.  At  higher  altitudes,  as  the  flow  becomes  more  rarefied,  the 
atmospheric  gases  can  penetrate  the  exhaust  plume  and  thereby  cause 
chemical  reactions  to  take  place  between  the  high  energy  gases  with  sub¬ 
sequent  emission  of  significant  additional  radiation.  The  need  for  infor¬ 
mation  on  the  basic  plume  structure  and  the  aerodynamics  of  this  mixing 
region  was  the  reason  for  the  present  experimental  study. 

Previous  theoretical  studies  have  been  made  (Refs.  2  through  7) 
of  rocket  exhaust  plumes  at  high  altitudes.  These  were  primarily  con¬ 
cerned  with  the  prediction  of  the  plume  size  and  shape.  Some  of  these 
(e.g. ,  Refs.  2  and  3)  were  semiempirical  in  that  simplified  solutions 
for  the  plume  boundary  were  obtained  by  fitting  equations  for  the  density 
variations  to  solutions  obtained  by  the  method  of  characteristics. 

Some  experimental  data  had  been  obtained  from  plumes  exhausting 
from  a  rocket  into  a  vacuum  or  into  a  static  pressure  region.  There 
was  almost  a  complete  absence  of  experimental  results  on  the  inter¬ 
action  of  a  plume  and  a  high  speed  rarefied  gas  stream.  The  limited 
data  that  had  been  obtained  comprised  photographic  and  spectrographic 
results  from  actual  rocket  flights  through  the  atmosphere. 

The  present  work  wag  carried  out  to  determine  the  density  varia¬ 
tion  of  the  exhaust  and  the  ambient  gases  across  the  plume,  as  well  as 
the  size  and  shape  of  the  free-stream  and  plume  interaction  region. 

The  effect  of  parameters  such  as  free-stream  Mach  number  and  density, 
exhaust  gas  properties,  and  rocket  nozzle  geometry  on  plume  char¬ 
acteristics  was  investigated. 


SECTION  II 

BASIC  THEORY  AND  SELECTION  OF  TEST  CONDITIONS 


2.1  GENERAL  FEATURES  OF  A  PLUME 

The  general  features  of  a  rocket  plume  at  high  altitude  are  shown 
in  Fig.  1.  The  plume  is  formed  by  the  ejqpansion  of  the  rocket  exhaust 
gases  into  the  rarefied  atmosphere.  The  oncoming  free-stream  gas 
causes  the  plume  gas  to  be  swept  downstream.  At  lower  altitudes  than 
that  for  which  Fig.  1  has  been  drawn,  the  free-stream  pressure  is 
higher,  ahd~the  plum'e“s‘iz‘e  consequently' smaller. 
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The  pLume  often  induces  flow  separation  on  the  vehicle  forward  of 
the  exit  plane  of  the  nozzle,  A  shodk  is  formed  in  front  of  the  separated 
region,  which  interacts  with  the  shock  from  the  nose  of  the  vehicle. 

The  prominent  bow  shock  is  initiated  by  this  interaction,  but  the  bow 
shock  is  influenced  primEiriLy  by  the  plxune  itself,  which  presents  the 
appearance  of  a  blunt  body  to  the  free  stream.  For  the  altitude  range 
of  interest  in  the  present  study,  the  plume  is  considerably  larger  than 
the  vehicle  so  that  the  details  of  the  interactions  of  the  vehicle  generated 
shocks  are  only  of  minor  importance. 

An  additional  shock  is  formed  within  the  plume,  because  of  the  coa¬ 
lescing  of  pressure  waves  from  the  curved  inviscid  contact  surface. 

This  annular  shock  tends  to  dose  on  itself  farther  downstream,  leading 
to  a  normal  shock,  called  the  Mach  disk,  at  some  point  in  the  flow. 
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Between  the  bow  shock  and  the  internal  shock  is  a  region  within 
which  the  plume  gas  and  the  free-stream  gas  are  in  contact.  It  is 
often  assumed,  for  simplicity,  that  the  two  gases  are  separated  by  an 
infinitesimally  thin  layer  called  the  inviscid  contact  surface.  Actually, 
diffusion  between  the  two  gases  must  take  place  so  that  there  is  a  mix¬ 
ing  region  of  finite  thickness  throughout  which  the  two  gases  are  pres¬ 
ent  in  varying  amounts. 

There  are  at  least  four  regions  of  the  flow  field  that  can  contribute 
to  the  total  plume  radiation  (Fig,  1): 

1.  The  region  of  hot  exhaust  gas  near  the  exit  plane  and 
farther  downstream. 

2.  The  region  of  exhaust  gas  that  has  passed  through  the 
internal  shock  or  Mach  disk  with  a  consequent  increase 
in  the  static  temperature. 

3.  The  atmospheric  gas  heated  by  the  strong  bow  shock. 

4.  The  mixing  region  between  the  two  gases  in  which 
chemical  reactions  between  the  components  of  the 
two  gases  can  take  place. 

The  purpose  of  the  experimental  study  reported  herein  was  to  es¬ 
tablish  the  fluid  mechanics  aspects  of  the  last  three  of  these  regions. 
This  was  done  by  determining  impact  pressure  distributions  through¬ 
out  the  plume  to  establish  the  relative  thicknesses  of  the  regions, 
and  by  determining  with  an  electron  beam  the  distributions  of  tempera¬ 
ture  and  density  of  the  gases,  with  the  emphasis  on  the  mixing  region. 

It  was  not  possible  to  use  an  actual  rocket  in  these  studies,  so  direct 
measurements  of  radiation  were  not  made.  Instead,  the  gases  were 
chosen  for  convenience  in  making  the  necessary  measurements. 

Thus,  nitrogen  was  used  for  the  free-stream  gas,  and  helium  for  the 
exhaust  gas,  because  the  electron  beam  techniques  for  determining 
densities  and  temperature  are  reasonably  well  established  for  these 
gases.  Carbon  dioxide  was  also  used  for  the  esdiaust  gas  for  some 
runs,  in  order  to  check  the  correctness  of  the  scaling  parameters  used, 
but  for  this  gas  determination  of  the  density  and  temperature  by  the 
electron  beam  technique  is  not  presently  practical. 


2.2  SIMPLIFIED  DESCRIPTIONS  OF  THE  FLQW_EIELD. 

The  flow  field  of  a  rocket  exhausting  into  a  vacuum  or  a  constant 
pressure  region  has  received  considerable  study.  The  case  of  a  rocket 
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exhausting  into  a  moving  free  stream  has  received  less  study.  Within 
the  assumption  of  inviscid  continuum  fluid  mechanics,  exact  solutions 
may  be  obtained  for  specific  cases  by  use  of  one  of  the  many  axisym- 
metric  characteristics  computer  codes  available.  There  is  a  place, 
however,  for  simpler  analytical  solutions.  Such  a  solution  has  been 
given  in  Ref.  2.  The  results  given  there  are  used  here  to  determine 
an  approximate  relation  for  the  location  of  the  inviscid  boundary  be¬ 
tween  the  plume  gas  and  the  free-stream  gas.  The  purpose  in  pre¬ 
senting  these  results  is  not  to  obtain  a  result  which  necessarily  gives 
the  proper  location,  but  to  show  what  the  basic  similarity  variables  are. 

Consider  first  the  case  of  a  rocket  e:daaustiag  into  a  vacuum  or  a 
constant  pressure  region,  as  shown  in  Fig.  2.  The  flow  field  within 
the  nozzle  is  normally  assumed  to  be  a  simple  source  flow,  or  a  uni¬ 
form  parallel  flow,  although  the  flow  field  within  an  actual  nozzle  may 
be  considerably  more  complicated.  At  the  corners  of  the  nozzle,  the 
flow  will  undergo  a  Prandtl-Meyer  expansion  through  an  an^e  suffic¬ 
ient  to  match  the  external  pressure.  The  bounding  streamline  will 
then  ciirve  inward,  because  of  the  two-dimensional  nature  of  the  flow. 
The  curvatime  of  this  line  leads  to  a  coalescing  of  pressure  waves 
and  the  generation  of  the  internal  shock. 


Fig.  2  Flow  Field  for  Expansion  into  a  Region  of  Uniform  Pressure 
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The  simplified  result  for  the  portion  of  the  plume  gas  within  the 
internal  shock  is  based  upon  the  observation  that  as  the  limiting  velocity 
is  reached,  Ug  approaches  zero.  It  then  follows  from  the  continuity 
equation  that 


pjUrr^  =  function  of  d 


It  must  be  emphasized  that  this  result  is  valid  only  for  Up  near  to 
^lim’  form  of  the  function  was  selected  on  the  basis  of  analytical 

simplicity,  and  by  fitting  a  number  of  exact  numerical  solutions.  The 
form  chosen  in  Ref.  2  and  used  here  is 

=  Kie-^^  (1  -  cos 

with  Ki  and  X  constants.  References  uses  the  function 

1  -  cos  ^im 

but  Eq.  (1)  is  said  to  provide  a  better  fit.  The  constants  Kj  and  X  are 
evaluated  by  requiring  that  the  integrated  forms  of  the  continuity  and 
momentum  equations  hold.  The  mass  flow  rate  through  the  surface  of 
constant  radius  r  is 


rii  ~  2?7  /  PiUfP*  sin  0  d9 
0 

Using  the  distribution  of  Eq.  (1)  gives 

m  =  (r7^''"Ki/A)  erf  (A(l  -  cos  ^nn,)) 

Practical  values  of  X  are  fairly  high  (> 5)  thus  the  value  of  the  error 
function  is  very  close  to  unity.  Therefore 

m  «  Ki/\  (2) 

This  must  be  equal  to  the  mass  flow  rate  of  plume  gas  through  the 
nozzle  throat 


m  =  p-U*A*  (3) 

Setting  the  error  function  equal  to  unity  is  necessary  for  reasons  other 
than  simplifying  the  resulting  equations.  The  flow  field  interior  to  the 
internal  shock  is  independent  of  the  boundary,  so  that  ©lim  should  not 
appear  in  the  final  result.  An  alternate  view  is  that  the  distribution  of 
mass  flow  rate  is  very  concentrated  near  the  centerline,  and  the  mass 
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flow  rate  in  the  outer  portions  of  the  plume  is  small.  Note  from  Eq.  (2) 
that  for  a  given  mass  flow  rate,  X  must  be  a  specific  constant.  Use  of 
a  variable  X  as  done  in  Ref.  2  is  not  justified. 

The  momentum  equation  is  integrated  along  the  exit  plane  and 
along  a  surface  in  the  flow  field  of  constant  r.  If  F  is  the  vacuum 
thrust  of  the  rocket,  then, 

Q 

F  =  Ur  cos  8  sin  6  dd 

0 

Letting  ®  and  setting  Uj,  to  gives 

F  =  KjU[5„  (1  -  i/v‘;rx) 

“  m  U  lim  U  *“  1,  \  ^7  X) 

The  thrust  coefficient,  Cp,  is  defined  by 

F  =  Cf  A*  Pp 


(4) 


(5) 


and  the  maximum  thrust  coefficient,  Cp^,  by 

-  Cf„  A*p^ 

=  mUiim 


From  these 


A  =  1.  vF(I 


For  an  ideal  gas 


^Fa,  y  _  1  ^y  + 


Y+  ^ 
2(y-ll 


(6) 


(7) 


(8) 


and 


IIZ 

~  y  2 


y-l  Med  +  l/yM.,-‘) 


/  y  -  I 


(9) 


Figure  3  gives  a  plot  of  X  versus  the- nozzle  area  ratio  for  given 
values  of  y. 
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Fig.  3  Plume  Spreading  Parameter,  X 


Now  the  problem  of  the  location  of  the  inviscid  contact  surface  is 
considered.  Figure  4  presents  a  simplified  sketch.  Basically,  the 
slope  of  the  contact  surface  at  a  location  r,  d  is  assumed  to  be  that  for 
which  the  pressure  on  each  side  of  the  contact  surface  is  the  same. 
Thus,  the  pressure  at  a,  which  is  determined  by  turning  the  free 
stream  through  the  angle  is  equal  to  the  pressure  at  2,  which  is 
determined  by  turning  the  plume  gas  through  the  angle  0  -  |3.  To 
-formulate  the  equations,  several  assumptions  are  made.  First,  the 
contact  surface  is  infinitesimally  thin,  and  the  regions  between  the  bow 
shock  and  the  internal  plume  shock  are  very  thin.  Second,  the  influence 
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Fig.  4  Simplified  Sketch  of  Plume  Contact  Surface 
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of  the  flow  over  the  body  on  the  plume  may  be  neglected.  Third,  the 
simple  Newtonian  pressure  law  may  be  used.  The  equality  of  the  pres¬ 
sures  then  leads  to 


P<»  +  =  Pi  +  />iUi^8in^(0  -  p) 


(10) 


The  next  assumption  is  that  the  distribution  of  Eq.  (1)  holds  within 
the  plume.  This  assumption  is  not  a  particularly  good  one,  because, 
as  mentioned  above,  the  distribution  given  by  Eq.  (1)  may  not  be  accu¬ 
rate  at  the  outer  edges  of  the  plume.  Nonetheless,  the  result  is  used, 
because  the  purpose  here  is  to  develop  the  basic  similarity  variables, 
not  determine  the  shape  of  the  plume  for  specific  cases.  Thus 

Pi^i  =  U,5n,  (11) 

r 


In  addition,  the  pressure  pj^  is  neglected.  By  using 

iIt 


sin  /3 


and 


sin(0  -  ft) 


r 


the  result  becomes 


Now  let 


(12) 
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and  the  result  becomes 


i  *  yfAj 

(1 


W2 


-4=  m 

2\ir 


The  boundary  conditions  to  be  applied  to  this  equation  are 


at  ^  W  = 


A  J  o  ^ 

r*  V  y^oo 


gpo 

Cf^Pc 


(13) 


(14) 


in  which  y^  is  the  exit  radius  of  the  nozzle.  Note  that  the  origin  of 
the  coordinate  system  has  been  located  on  the  centerline  at  the  inter¬ 
section  of  the  extension  of  the  slope  of  the  contact  surface  at  the  cor¬ 
ner.  This  has  been  done  to  make  the  equations  somewhat  simpler j 
although  in  general  the  distance  y^  is  small  in  comparison  with  the 
pliime  size. 


The  maximum  turning  angle  is  established  by  a  pressure  balance 
at  the  corner.  A  sketch  is  shown  in  Fig.  4.  The  turning  angle  of  the 
nozzle  flow  is  pressxjre  is  then  obtained  from  the  isen- 

tropic  taldes,  i.  e. , 

p  s  —  -  'K  )  Pe 

Po 


in  which  v  is  the  PrandtL-Meyer  function.  Setting  this  pressure  equal 
to  the  pressure  attributable  to  the  free  stream  (and  neglecting  p^^  gives 

~  pr  —  bKMr)  +  fJm  -  "v) 


or 


^  -  0^) 
2j3 _ 


2  «in 


Pc 


(15) 


For  a  given  nozzle,  and  given  7,  the  value  of  therefore  depend¬ 

ent  upon  the  ratio  of  the  chamber  pressure  to  the  free-stream  dynamic 
pressure  (the  result  can  also  be  used  for  the  case  of  no  flow  by  setting 
the  pressure  equal  to  p^. 


Combining  Eqs.  (12),  (13),  and  (14)  leads  to 


y-y. 


yMne 

1  +  nin“  6 


Cp  p, 

^  m  e 


il 
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in  which  0^^  is  a  function  of  Pq/ Usually  the  influence  of  is 
minor  and  may  be  neglected.  Also,  for  simplicity,  yg  may  be  neg¬ 
lected.  Then 

^  m 


This  yields  the  basic  similarity  parameters.  The  parameters  q^,,  pg, 
and  r*  can  be  absorbed  into  a  y  scaling  parameter,  defined  by 


L 


(17) 


2 

The  primary  remaining  variables  are  X  and  7M„  .  Of  these,  X  is  of 
primary  importance  and  of  secondary  importance.  In  fact,  if 


yM^  sin ”5  >  >  1 


then 


y  -y. 


d 


(18) 


The  parameter  X  shoiild  not  be  absorbed  into  the  y  scaling  in  this  case, 
because  the  parameter  X  also  appears  in  the  function  f(e),  Eq.  (1).  The 
particular  form  of  Eq.  (1)  is  such  that  an  analytical  integration  does  not 
appear  possible.  However,  if  the  approximation 

2 

1  —  cos  0  1  —  $ 

V 


is  made  (see  Fig.  5),  then  the  result  may  be  integrated  to  give 

(^) "  ^  «  -  ""  ®">)  -  (^r 


or  by  setting  the  value  of  the  first  error  function  to  unity,  as  done 
earlier,  gives 


Pc  V  '  / 


(19) 


The  value  of  Eq.  (19)  is  not  in  the  accuracy  with  which  it  represents 
more  exact  solutions,  but  it  is  the  explicit  form'in  which  the' basic 
scaling  parameters  appear.  From  this  it  can  be  seen  that  the  main 
parameter  is  X.  Other  parameters  appear  in  groupings  such  that  it- 
is  not  necessary  to  match  actual  values,  and  ratios  of  parameters  may 
be  matched. 
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The  scaling  parameters  considered  up  to  this  point  have  been 
established  on  the  basis  of  inviscid  continuum  fluid  mechanics.  Some 
other  factors  of  importance  are  associated  with  viscous  effects  and 
low  density  effects.  These  are  important  primarily  in  the  mixing 
region.  The  simulation  parameters  of  this  region  are  considerably 
more  complicated,  and  no  attempt  will  be  made  here  to  establish  the 
specific  parameters.  A  preliminary  look  was  obtained  by  testing  at 
two  different  density  levels  by  using  two  wind  tunnels. 


2.3  SELECTION  OF  TEST  CONDITIONS 

The  tests  were  riui  in  two  VKF  wind  tunnels:  Low  Density  Hyper¬ 
sonic  Wind  Tunnel  (M)  and  Supersonic  Wind  Tunnel  (D).  Descriptions 
of  these  tunnels  are  given  in  Sectidn  m.  The  operating  parameters  of 
these  tunnels  established  the  test  section  conditions,  as  given  in 
Table  I,  and  the  size  of  the  plume  simulator. 

For  a  number  of  reasons  it  was  desirable  to  use  helium  as  the 
plume  gas.  The  first  of  these  was  because  of  the  suitability  of  helium 
for  density  measurements  by  use  of  the  electron  beam  technique.  The 
second  was  that  helium  has  a  very  low  condensation  temperature  and, 
therefore,  is  less  likely  to  undergo  condensation  within  the  plume. 
Therefore,  the  helium  nozzle  was  scaled  to  give  the  same  value  of  X  as 
a  real  nozzle,  for  which  a  value  of  7  of  1.27  was  selected, 

For  the  full-scale  flight  nozzle,  the  following  was  assumed: 

Ae/A’'  =  25 
y  =  1.27 


Therefore, 


A  =  4.7 

For  a  helium  exhaust  from  a  model  nozzle, 

y  =  1.67 


and  for  the  same  value  of  A, 


Ae/.A*  =  1.61 

The  maximum  feasible  chamber  pressure  for  Tunnel  M  was  found  to  be 
5  psi  for  the  model  tested.  For  the  same  value  of  the  correspond¬ 
ing  chamber  pressure  for  flight  is  440  psi. 
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TABLE  I 

SUMMARY  OF  TEST  CONDITIONS 


I 


I 


Coiidit  Lon 


Tonne] 

Free-  stream 
Gas 


Free-stream 

Conditions 


J.I. 


P—  pa  I 
Psi 

Re*,  per  in» 


Altitude,  Km 


Kazzle 

Cab 

t 

Ag/A* 

Me 

d*.  In. 

'dg,  m. 

Pc-  Ps* 

Sm 

''  *  r+-\/cF,„Pc/q- 
Purpose 


‘Data 

Impact 
'Pressure 
Electrcir  Beam 
Photograpliy 
P,  T 


Pirsi  Condition 


JC5 

jAB 

k1 

kl 

Air 

N2 

N, 

— 

... _ 

..  .Second  Condition 

IIB 

Piiklu 

i£2_ 

lAC 

LftA 

JB 

U 

M 

M 

M 

M 

N2 

Air 

N2 

N2 

Nj 

N2 

18.00 
1.5  X  10'< 
3.  4  X  10'® 
200 
0.  12 
0.  024 
80 


1.27 

25 

4.  26 
12.9 
64.  7 
440 
41 
4.  7 
87  fl 


18.  15 

18.  IS 

I.  2  X  lO"'* 

1,  2  X  10'^ 

2.  7  X  10-2 

2.  7  X  10'^ 

1250 

1250 

0.  022 

0.  022 

1.  69 

Ie8l 

He 

COo 

1.67 

1.27 

1.61 

25 

2.  07 

4.26 

0.  315 

0.  03U 

0.  40 

0.  15 

5 

400 

42 

43 

4.5 

4.  7 

2.7  m. 

2.58  1 

5,0 

ID.  00 

1.5  X  10'® 

I,  5  X  10-* 

2.0  X  10'* 

1.  0  X  10-2 

2.9  X  lu8 

125 

0.  00021 

0.  12 

0.015 

0.039 

80 

He 

1.87 

1.2? 

1.61 

25 

2.  07 

4.  26 

0.  315 

12.9 

0.  40 

64.  7 

56 

53 

42 

36 

4.  5 

4.  7 

2.  94  in. 

54  ft 

>8.15  18,15 

1.2  X  10-'>  1.2  X  lo-l 

2.  7  X  10-2  2.7  X  10-2 

1250  1250 

0,  022  0.  022 

2.68  1.3fi 


18.  15 

18.15 

1.2  X  10-< 

1.2  X  10“^ 

2.  7  X  10-2 

2.7  X  10-2 

1250 

1250 

0.  022 

0.  022 

1.  88 

5.  39 

He 

1.67 

CO, 

1.27 

COj 

1.27 

He 

1.67 

1.  61 

25 

25 

25 

2.  07 

4.  26 

4.  26 

7.  09 

0.315 

0.  080 

0.  080 

0.  030 

0.  40 

0.  40 

0.40 

0.40 

2 

100 

50 

450 

37 

36 

32 

22 

4.5 

4.  7 

4.  7 

33.5 

1,  70  in. 

3.  44  in. 

2.  43  in. 

2.  31  in, 

Basic 

data 


Check 

SIIIIU- 

lalhMI 


Check 

effvct 


Basic 

data 


Check 

simu¬ 

lation 


Check 

Pe 

effect 


Check 

T 

effect 


N 


\ 


x  .X 

X  X 


X  X 

X  X 


\ 
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The  conditions  are  summarized  in  Table  I.  The  symbolism  for 
the  conditions  unfortunately  underwent  considerable  alteration  as  modi¬ 
fications  of  the  original  test  plan  were  found  to  be  necessary,  so  that 
the  final  designations  do  not  have  the  consistency  that  would  be  pres¬ 
ent  if  the  designations  had  been  made  later.  The  designations  have 
been  preserved  in  the  form  used  during  the  experimental  portion  of 
the  program. 

The  first  four  columns  list  conditions  for  tests  corresponding  to 
one  flight  condition.  Condition  ICS  was  the  corresponding  case  in 
Tunnel  M,  using  helium  as  the  plvime  gas.  Condition  lAB  was  run  in 
Tunnel  M  with  CO2  as  the  plume  gas,  in  order  to  check  the  basic  scal¬ 
ing  parameters.  Condition  HB  was  run  in  Tunnel  D  to  check  the  effect 
of  greater  free-stream  density. 

The  next  three  columns  correspond  to  the  first  three,  but  with  a 
reduced  chamber  pressure.  The  last  two  columns  represent  some 
extra  conditions. 

It  was  sometimes  impossible  to  achieve  the  desired  conditions  in 
the  tunnel.  The  reasons  were  various,  but  typically  were  the  inability 
to  pump  the  desired  mass  flow  rate  or  the  inability  to  achieve  the 
desired  maximum  pressure.  Therefore,  there  was  not  always  obtained 
in  practice  the  exact  correspondence  desired. 

The  various  conditions  are  discussed  in  order  in  the  following  para¬ 
graphs.  The  first  condition  (see  Table  I)  is  for  flight  at  an  altitude  of 
80  km  at  a  free-stream  Mach  number  of  18,  with  the  latter  chosen  to 
correspond  to  the  test  section  Mach  number  of  Tunnel  M.  The  nozzle 
geometry  was  given  earlier.  The  chamber  pressure  of  440  psi  was 
selected  to  give  the  same  value  of  obtained  in  condition  IC5 . 

Condition  ICS  was  run  in  Tunnel  M  with  the  helium  nozzle.  The 
chamber  pressure  was  set  at  the  maximum  value  that  could  be  main¬ 
tained  without  building  up  excessive  backpressure  in  the  system.  The 
size  of  the  nozzle  was  sd.ected  so  that  the  maximum  plume  diameter 
would  be  about  8  in. 

Condition  lAB  was  run  in  Tunnel  M  with  the  small  CO2  nozzle. 

The  intention  was  to  set  the  chamber  pressure  to  the  same  as  the  flight 
value  (440  psi),  but  limitations  on  the  system  made  it  necessary  to  run 
at  a  slightly  io^r” value.  '  '  ”  ■  ■  ■  “ 

Condition  IIB  was  run  in  Tunnd  D  with  the  helium  nozzle.  The 
primary  purpose  for  this  test  was  to  determine  the  rarefied  flow  effects. 
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The  free-stream  dynamic  pressure  in  Tunnel  D  was  approximately 
10  times  that  in  Tunnel  M,  and  the  nozzle  chamber  pressure  was 
therefore  also  set  at  about  10  times  the  value  for  M.  The  mean  free 
path  in  Tunnel  D  was  0.  01  times  that  in  Tunnel  M,  so  that  any  rarefied 
flow  effects  would  be  considerably  reduced. 

The  second  condition  was  selected  the  same  as  ICS,  but  with  a 
reduced  chamber  pressure  of  2  psi,  rather  than  5  psi.  Both  the 
helium  nozzle  (condition  IC2)  and  the  large  CO2  nozzle  (condition  lAC) 
were  used.  The  corresponding  flight  condition  was  a  free-stream 
Mach  number  of  10  at  80  km. 

Condition  lAA  was  run  with  the  CO2  nozzle  at  a  chamber  pressure 
of  50  psi.  Therefore,  conditions  lAA  (50  psi),  lAB  (400  psi),  and  lAC 
(100  psi)  allow  direct  comparisons  to  be  made  on  the  effect  of  pressure. 
Condition  IB  was  run  to  get  some  direct  information  on  the  effect  of  7. 
The  condition  was  run  with  the  same  nozzle  and  chamber  pressure  as 
condition  lAB,  but  with  helium  rather  than  COg  as  the  plume  gas. 

Characteristic  solutions  for  some  of  these  conditions  are  given  in 
Figs.  6a  through  d. 


SECTION  111 

DESCRIPTION  OF  EQUIPMENT 


3.1  TUNNEL  M 

The  experimental  work  was  performed  in  Tunnels  M  and  D  of  the 
von  Karman  Gas  Dynamics  Facility  (VKF).  Tunnel  M,  shown  pictorially 
in  Fig.  7a  and  schematically  in  Fig.  7b,  is  a  continuous,  arc-heated, 
low-density,  hypersonic  wind  tunnel  normally  using  nitrogen  as  the  test 
gas.  Pumping  is  provided  by  three  stages  of  air  ejectors  in  series 
which  exhaust  into  the  VKF  main  compressor  system  through  the  VKF 
Hypersonic  Wind  Tunnel  (C)  test  section.  This  arrangement  permits 
simultaneous  operation  of  these  two  tunnels,  or  either  can  be  operated 
alone.  Tunnel  M  consists  basically  of  the  following  major  components, 
in  streamwise  order: 
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1.  Rotating- arc -type  d-c  arc  heater  of  VKF  design  with 
a  power  supply  rated  at  200  kw  for  continuous  opera¬ 
tion.  Gas  is  injected  into  the  arc  heater  tangentially. 

2.  Cylindrical  settling  chamber  of  1.  5 -in.  diameter  and 
3 -in.  length. 

3.  Axisymmetric,  contoured  aerodynamic  nozzle  with 

0. 181-in.  throat  diameter  and  35. 1-in.  exit  diameter. 

A  6-  to  10-in. -diam,  gradient -free  test  core  extends 
from  25  in.  upstream  of  the  nozzle  exit  to  25  in. 
downstream  of  the  exit.  Table  II  gives  the  nominal 
Tunnel  M  operating  conditions  with  the  contoured 
nozzle. 

4.  Stationary  bulkhead  of  94-in.  diameter  which  supports 
the  nozzle,  probe  drive  and  support  unit,  pressure 
measuring  system,  and  external  force  balance  or 
model  support  base.  The  bxilkhead  contains  eight 
12-in.-diam  ports. 

5.  Cylindrical  8-ft-diam  test  chamber  which  moves  down¬ 
stream  to  allow  access  to  the  test  section,  models, 
and  probes. 

6.  Axisymmetric  diffuser  with  convergent  entrance,  con¬ 
stant  area  throat,  and  divergent  outlet.  Interchange¬ 
able  designs  are  available  for  different  test  configura¬ 
tions. 

7.  Downstream  heat  exchanger. 

8.  First  air  ejector  stage. 

9.  Isolation  valve. 

10.  Second  and  third  air  ejector  stages. 

11.  Connection  to  the  VKF  Tunnel  C  test  section. 

Backside  water  cooling  is  used  to  protect  the  arc  heater,  settling  cham¬ 
ber,  nozzle,  and  diffuser. 
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a.  Condition  IC5 
Fig.  6  Characteristics  Solutions 


c.  Condition  IC2 
Fig.  6  Continued 
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a.  Photograph 


C  TffT  UCTKM  • 


CLEVUION  VIEW  OF  TUNNEL  M 

b.  Elevation  View/ 
Fig.  7  Tunnel  M 
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TABLE  II 

NOMINAL  TUNNEL  M  NOZZLE  I8N2  OPERATING  CONDITIONS 


a.  Customary  Units 


m,  Ibj^/sec 
p  .  atm 
»K 
ho* 

pQ.  mm  Hg 

^  1 
Re„,  in.  "  ^ 

P«* 

T„,  “K 
U„,  ft/  sec 
Ibj^/ft^ 
in. 

qa>.  Ibf/ft2 
Re2j  in. 
Reo,  in. 


0. 0515 
19.  0 
2900 
1500 

2.  55 
18.  15 
1250 
6.  00 
45.  0 
8165 

3.  73  X  10“6 
0.  022 
3.90 

45.  0 
44.  0 


b.  International  System  of  Units 


in,  kg/ sec 

0.0233 

Pq.  N/m2 

1. 925  X  10 

To*  "K 

2900 

ho*  J/kg_ 

3.49  X  10® 

Po*  N/m^ 

340 

18.  15 

Re„,  m"^ 

4.  92  X  10^ 

p„,  N/m2 

0.  800 

T„,  “K 

45.  0 

U„,  m/sec 

2490 

P„,  kg/m^ 

5. 97  X  10"' 

5.  59  X  10"' 

,, 

q«>  N/m“ 

187 

Re2i  m“^ 

Rcq,  m”^ 

1770 

1730 
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3.2  TUNNEL  D 

Tunnel  D  is  an  intermittent,  variable  density  wind  tunnel  with  a 
manually  adjusted,  flexible-plate-type  nozzle  and  a  12-  by  12 -in.  test 
section.  A  picture  is  given  in  Fig.  8a  and  a  schematic  in  Fig.  8b. 

The  tunnel  operates  at  Mach  numbers  from  1.  5  to  5  at  stagnation  tem¬ 
peratures  up  to  about  540®R. 


3.3  PLUME  GENERATION  SYSTEM 

Three  rocket  simulators,  shown  in  Fig.  9,  were  used  to  generate 
the  various  plumes  studied  in  Tunnel  M.  One  of  these  (d*  =  0.315  in. ) 
was  used  in  Tunnel  D.  The  0.  315-in.-throat-diam  model  had  a  nozzle 
expansion  ratio  of  1.  6,  whereas  the  0.  030-  and  0.  080-in.-throat-diam 
models  had  expansion  ratios  of  25.  All  had  a  10-deg  half-angle  coni¬ 
cal  expansion.  Dimensions,  given  in  Table  HI,  were  kept  as  small 
as  practicable  to  minimize  the  effect  of  the  simulator  itself  on  the  free 
stream.  The  chamber  diameter  of  the  0.  315-in.-throat-diam  model 
is  of  the  same  order  as  the  throat  diameter.  Thus  it  was  machined 
with  a  perforated  conical  baffle  to  promote  flow  uniformity  entering  the 
nozzle.  A  water-cooled  strut,  which  supported  the  simulator  in  the 
test  section,  and  gas  supply  and  pressure  lines  extended  up  from  the 
model  support  base  and  attached  to  the  lower  wall  of  the  simulator 
chamber.  A  typical  installation  in  Tunnel  M  with  impact  pressure 
probes  is  shown  in  Fig,  10a. 

The  plume-gas  supply  system  is  shown  schematically  in  Fig.  11. 
Helium  (99.  995-percent  purity)  was  supplied  from  a  41,000-ft3 
(standard  temperature  and  pressure,  STP)  transport  trailer  initially 
at  2000  psi.  The  small  quantity  of  argon  used  (99.  99  +  -percent  purity) 
came  from  2000  psi,  220-ft3  (STP)  cylinders.  Carbon  dioxide 
(99.  5-percent  minimiim  purity)  was  obtained  as  flash-off  from  a  mani¬ 
folded  bank  of  nine  50 -lb  liquid  CO2  cylinders  at  700  psi.  A  manually 
operated  Hoke  pressure  regulator  in  series  with  two  manual  throttling 
valves  was  used  to  regulate  gas  flow  to  the  heater  and  maintain  the 
desired  pressure  in  the  rocket  simulator  chamber. 

In  Tunnel  M  the  plume  gas  was  heated  to  approximately  two -thirds 
the  desired  chamber  temperature  by  passing  it  through  two  parallel- 
connected,  single-pass,  hot-core,  Kanthal®-element  resistance 
heaters  installed  for  these  experiments.  Final  heating  was  provided 
by  the  tunnel  flow  (T^  =  2900*10  impinging  on  the  simulator  nose  and 
supply  tube.  Desired  chamber  temperaturje  was  maintained  by-  regu¬ 
lating  the  heater  output  by  means  of  a  manually  operated  control 
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b.  Schematic 
Fig.  8  Tunnel  D 
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Fig.  9  Photograph  of  Plume  Generators 
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b.  Tunnel  D 

Fig.  10  Simulator  and  Impact  Pressure  Probe  Installation 
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TABLE  III 

PLUME  SIMULATOR  DIMENSIONS 


! 

Throat 

Nozzle 

Chamber 

Diametef, 

Diameter  j 

Diameter, 

in. 

d]^,  in. 

d2>  in. 

0.  315 

0.  400 

0.406 

0.  080 

0.  400 

0.  375 

0.  030 

1  0.  150 

0.  150 

L>2>  L,3, 

in.  in. 

0.  638  0.  859 

0.284  0.687 

0. 284  0. 343 


Thermocouple 


Heated  Plume 
Gas  In 


To  Instrument  for 
Measurement  of  T, 


To  Pressure  Transducer 
“for“M'ea'snr  ement“of-p^ - 
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rheostat.  For  Tunnel  D,  no  heater  was  used.  Figure  10b  shows  the 
simulator  and  impact  pressure  probes  installed  in  Tunnel  D. 

Estimates  of  the  boundary-layer  thickness  at  the  exit  plane  of  the 
simulator  nozzles  give  values  of  about  0.  030  in.  for  conditions  IC2 
and  lAC,  The  corresponding  value  of  6*  is  about  0.  020  in.  These 
values  are  felt  to  be  sufficiently  low  that  there  should  be  no  signifi¬ 
cant  viscous  effects. 


3.4  TUNNEL  CONDITIONS  AND  PRESSURE  MEASURING  EQUIPMENT 

The  instrumentation  apparatus  is  shown  in  Figs.  12a  and  b.  The 
tunnel  free- stream  conditions  were  monitored  from  measurements  of 
reservoir  pressure  and  test  section  impact  pressure.  The  third 
measurement  was  the  reservoir  temperature  in  Tunnel  D,  which  was 
measured  with  a  thermocouple,  and  the  mass  flow  rate  in  Tunnel  M, 
from  which  a  reservoir  temperature  was  determined. 

In  Tunnel  M  impact  pressures  in  the  free  stream  and  in  the  plume 
were  meastired  with  water-cooled  probes  connected,  respectively,  to 
3-  and  30-mm  MKS  Baratron®  pressure  meters.  The  flat-faced  free- 
stream  probe,  located  approximately  even  with  the  nose  of  the  rocket 
simulator  and  3  in.  below  the  nozzle  centerline,  was  0.  2 5 -in.  OD  and 
had  a  0.  093-in.-ID  pressure  tube.  The  externally  chamfered  plume 
probe,  with  a  0.  121-in,-diam  nose  and  a  0.  063-in.-diam  orifice,  was 
mounted  on  the  probe  carrier  so  that  it  could  be  traversed  axially  and 
radially  through  the  plume  flow  field.  Two  10 -turn  potentiometers 
were  used  to  indicate  probe  position.  The  Baratron  pressure  meters 
were  calibrated  each  day,  using  a  precision  oil  micromanometer  and 
a  McLeod  gage,  under  the  same  ambient  conditions  which  they  exper¬ 
ienced  during  the  test.  All  other  pressure  instruments  were  cali¬ 
brated  before  the  test. 

Digital  voltmeters  plus  various  analog  readout  devices  were  used 
to  monitor  various  parameters  during  the  test.  All  analog  data  were 
fed  into  an  analog-to-digital  converter  and  then  recorded  on  paper  tape 
fob  subsequent  computer  processing. 

In  Tunnel  D,  two  probes  were  attached  to  a  common  positioning 
mechanism  (Fig.  13).  One  of  these  probes  was  angled  to  bring  the 
probe  in  closer  alignment  with  the  plume  flow.  There  was  some  over¬ 
lap  in  the  coverage  of  the  plume  by  these  two  probes. 
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b.  Tunnel  D 
Fig.  *12  Concluded 
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Fig.  13  Tunnel  D  Simulator  and  Impact  Pressure  Probes  during  Typical  Run 
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The  impact  pressures  in  Tunnel  D  were  measured  on  a  Consolidated 
Electrodynamics  Corporation  (CEC)  electromanometer  pressure  system 
using  15-psid  precision  pressure  balances  referenced  to  a  vacuum. 

The  rocket  simulator  nozzle  chamber  pressure  was  measured  on 
the  same  system  using  a  60-psid  precision  pressure  balance  referenced 
to  a  vacuum. 


3.5  ELECTRON  BEAM  APPARATUS 

Measurement  of  the  rotational  and  vibrational  temperatures  of 
nitrogen  and  number  density  of  both  nitrogen  and  helium  throughout  the 
plume/free  stream  region  was  accomplished  in  Tunnel  M  using  the 
electron  beam  fluorescence  technique  described  in  Appendix  II.  The 
apparatus  required  for  injecting  the  high  energy  beam  of  electrons  into 
the  low  density,  hypersonic  test  gas  is  simple  in  nature,  requiring  only 
an  electron  source,  power  supply,  and  differential  pumping  system.  A 
photograph  of  the  electron  beam  apparatus  installed  in  Tunnel  M  is 
shown  in  Fig.  14.  The  source  of  electrons  was  a  GE  Type  33  television- 
type  electron  gun,  with  an  oxide-coated  cathode,  modified  for  50-kv  opera¬ 
tion.  The  accelerator  grid  had  an  enlarged  opening  of  approximately 
3/16-ln.  diameter,  allowing  beam  currents  of  up  to  2  ma  to  be  injected 
into  the  flow.  The  power  supply  for  the  filaments  consisted  of  batteries, 
yielding  6  to  18  vdc.  To  obtain  the  required  energy  of  approximately 
20  kev,  the  anode  of  the  gun  was  grounded  and  the  cathode  connected 
directly  to  the  negative  output  of  a  0-  to  50-kv,  0-  to5-ma  power  supply. 
Both  the  grid  and  the  preaccelerator  bias  voltages  were  supplied  by 
batteries.  Figure  15  shows  the  circuit  diagram  of  the  gxm  power  supply. 
The  gun  Itself  was  mounted  within  the  tunnel  in  a  leaktight  insixlated 
housing  which  was  maintained  at  1-atm  pressure  of  air,  as  were  the 
tubes  leading  the  power  supply  cables  into  the  tunnel;  this  arrangement 
prevented  arcing  to  the  tunnel  wall.  The  beam  was  magnetically  focused 
and  deflected  to  coincide  with  an  orifice  of  approximately  0.  050 -in. 
diameter  opening  into  the  test  section  containing  the  flow.  To  maintain 
the  pressure  within  the  electron  gun  at  approximately  10“^  torr,  as 
required  for  stable  and  long-life  operation,  a  differential  pumping  sys¬ 
tem  was  employed.  Referring  to  Fig.  16a,  it  can  be  seen  that  the  elec¬ 
tron  gun  was  connected  to  the  flow  by  a  24-in. -long  stainless  steel  drift 
tube  of  2 -in.  diameter;  connected  to  this  drift  tube  by  means  of  a  flex¬ 
ible  bellows  tube  was  a  pumping  line  leading  to  a  4 -in.  diffusion  pump 
located  outside  the  tunnel  (Fig.  16b).  A  solenoid-operated,  2 -in.  gate 
valve  was  located  at  the  flow  field  end  of  the  drift  tube. 
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Fig.  14  Photograph  of  Electron  Beam  Apparatus  in  Tunnel  M 
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Fig.  15  Electron  Gun  Power  Supply  Schematic 
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Rg.  16  Electron  Beam  Schematics 
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b.  Electron  Beam  Vacuum  Systems 
Fig.  16  Concluded 
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Measurement  of  the  beam  cxarrent  injected  into  the  flow  and  the 
associated  difficulties  are  described  in  Appendix  n.  For  vibrational 
and  rotational  temperature  measurements  a  precise  value  of  the  current 
is  unimportant  except  for  possible  secondary  effects  which  cannot  at 
this  time  be  included  in  the  equations  given  in  Appendix  H.  For  this 
reason  the  value  of  cimrent  used  for  all  temperature  data  was  a  nomi¬ 
nal  value  as  registered  by  an  aluminxim  flat-plate  collector  of  approxi¬ 
mate  dimensions  12  by  18  in.  located  on  the  tunnel  floor  approximately 
4  ft  from  the  orifice.  The  accuracy  of  such  a  measurement  is  indeed 
questionable  because  of  both  boundary  layer  and  secondary  ejected 
electron  effects;  however,  the  value  of  current  was  constantly  moni¬ 
tored  to  ensure  minimum  current  fluctuations  and  drift,  thereby  de¬ 
creasing  the  imprecision  of  the  temperature  measurements.  As  can 
be  seen  from  Eq.  (11-2),  knowledge  of  the  value  of  the  beam  current 
is  critical  for  an  accurate  density  measurement;  for  that  reason  it  was 
deemed  essential  that  an  additional  current  determination  be  made, 
specifically  at  the  point  of  injection  of  the  beam  into  the  flow,  thereby 
minimizing  the  errors  attributable  to  beam  spreading  caused  by  elastic 
and  inelastic  scattering.  To  do  this,  a  collector,  to  be  described  shortly, 
was  moimted  onto  a  traversing  mechanism  powered  by  a  24-vdc  power 
supply.  Before  obtaining  a  density  datum  point,  the  auxiliary  collector 
was  traversed  axially  imtil  it  was  located  directly  beneath  the  gun  orifice 
(ETg.  16a),  The  current  at  the  injection  point  was  found,  and  the 
auxiliary  collector  was  withdrawn  to  allow  the  traversal  of  the  flow  by 
the  beam.  Such  a  collector  must  meet  one  primary  requirement,  namely 
that  it  measure  all  current  incident  upon  it.  Shown  in  Fig.  17  is  the 
geometrical  configuration  of  the  collector  which  was  constructed  of  a 
2 -in.  copper  elbow  with  a  reflecting  plane  soldered  in  place  to  give  an 
effective  increase  in  the  length- to -diameter  ratio,  thereby  minimizing 
the  loss  of  ejected  secondary  electrons.  The  location  of  this  current 
measurement  site  was  approximately  24  in.  above  the  tunnel  center- 
line.  The  individual  values  of  current  from  eadb  of  the  two  collectors 
were  displayed  and  read  from  a  digital  voltmeter. 

The  optical  radiation  produced  as  a  result  of  electron  beam  excita¬ 
tion  was  gathered  by  a  48 -in. -focal -length  lens  which  has  an  approximate 
f/number  of  12  and  was  focused  onto  the  slit  of  a  Jarrell -Ash  1-meter 
scanAing  spectrometer  of  approximate  f/8.  6  as  is  shown  in  Figs.  18 
and  19.  The  102-  by  l02-mm  diffraction  grating  used  had  1180 
grooves/mm  and  was  blazed  at  5000  A  ,  For  the  measurement  of  all 
qu antiti es “the '’obj'ectli'istance  and'fET^imag^  distance  tb“tRe"Iens  were 
approximately  48  and  24  in. ,  respectively,  yielding  a  lens  magnifica¬ 
tion  factor  of  1/2.  Temperattire  measurements  were  performed  with 
the  spectrometer  slit  aligned  parallel  to  the  beam,  and  a  dove  prism  ' 
was  employed  for  the  density  measirnements  to  rotate  the  beam. 
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Fig.  17  Photograph  of  Current  Collector 


Fig.  18  Photograph  of  Spectrometer 
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perpendicular  to  the  slit.  Detection  was  accomplished  using  an 
EMI  6256B  photomultiplier  tube  with  the  output  processed  hy  a 
PAR  HR-8  lock-in  amplifier.  The  optical  signal  to  the  spectrometer 
was  chopped  by  a  PAR  Model  125  chopping  wheel  at  340  Hz,  and  the 
reference  signal  to  the  lock-in  amplifier  was  supplied  by  the  chopping 
wheel  assembly.  Rotational  temperature  data  were  obtained  b;  re¬ 
cording  the  lodk-in  output  on  a  stripchart  recorder.  For  both  vibra¬ 
tional  temperature  and  density  measurements  the  output  of  the  lock-in 
was  fed  to  an  active  electronic  integrator,  the  schematic  of  which  is 
shown  in  Pig.  20,  which  had  a  time  constant  of  1  sec  and  an  accuracy 
of  ±1  percent  for  integration  times  less  than  or  equal  to  200  sec.  The 
output  of  the  integrator  was  then  displayed  on  a  digital  voltmeter  from 
which  the  final  data  were  obtained.  Necessary  linearity  checks  of  the 
system  had  been  performed  before  the  test.  j 

To  achieve  the  axial  scanning  capability  required  for  high  quality 
photography  of  the  flow  field  as  well  as  fine  adjustment  to  the  beam 
position,  the  entire  electron  gun  housing  was  mounted  on  a  probe 
carrier  driven  by  a  spline  drive.  The  axial  position  of  the  beam  was 
determined  by  displaying  a  potentiometric  output  on  a  digital  volt¬ 
meter  calibrated  to  give  3000  counts  per  10  in.  Rotation  in  the  plane 
perpendicular  to  the  axial  direction  was  also  possible  and  was  used 
when  aligning  the  beam.  Radial  scanning  was  performed  by  mounting 
the  lens  and  spectrometer  on  a  ’’welder’s -table”  which  was  capable  of 
being  raised  or  lowered  by  crank,  thereby  moving  the  spectrometer 
and  optics  system  along  the  radial  direction  of  the  beam.  To  ensure 
that  the  electron  beam  {to' first  approximation)  intersected  the  tunnel 
centerline,  an  object  was  located  on  the  centerline  and  was  focused 
into  view  using  a  levelling  telescope  on  the  side  of  the  tunnel  opposite 
the  spectrometer  (see  Fig.  18).  A  plumb  bob  was  then  suspended 
from  the  electron  beam  orifice  and  the  beam  mount  was  traversed 
until  the  plumb  bob  intersected  the  centerline,  which  was  also  indi¬ 
cated  by  the  image  in  the  level.  The  lens  was  then  adjusted  mitil  the 
entrance  slit  of  the  sp)ectrometer  came  into  focus  with  the  level,  and 
at  that  point  it  was  certain  that  (1)  the  electron  beam  intersected  the 
centerline  and  (2)  the  spectrometer  entrance  slit  was  focused  onto  the 
beam.  After  the  magnetic  field  coil  was  turned  on  during  tunnel  opera¬ 
tion  a  final  adjustment  of  the  axial  and  angular  position  of  the  beam  was 
required,  but  the  change  was  rather  small  and  did  not  require  refocusing. 
The  radial  position  of  observation  was  determined  by  mounting  a  scale 
-on  the  spectrometer-tabLe,'  using-the-level-,—a-centerline -reference  posi= - 
tion  was  determined  on  the  scale,  which  had  a  1-mm  least  count,  and 
the  table  was  raised  until,  as  observed  through  the  level,  a  change  of 
position  of  predetermined  magnitude  had  occurred.  * 
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3.6  PHOTOGRAPHIC  APPARATUS 

For  flow  visualization  photographs,  a  camera  was  mounted  on  the 
surveyor's  stand  that  normally  held  the  levelling  telescope  used  for 
electron  be  am -spectrometer  alignment.  The  camera  was  a  Super 
Speed  Graphic®  with  a  135-mm  Graflex  Optar®f/4,  7  lens  and  a 
Graphex-Wollensak®  shutter  assembly. 

For  determination  of  proper  film  exposure  time,  Polaroid  Type  52 
(black  and  white)  and  Polaroid  Type  108  (color)  were  normally  used. 

The  film  used  for  data  point  photographs  was  either  Ektacolor®Type  Lf 
or  High  Speed  Ektachrome®  B  (processed  as  a  negative). 

Filters  were  frequently  used  with  data  point  photographs  to  prevent 
the  predominant  radiation  of  some  molecular  species  excited  by  the 
electron  beam  from  exposing  the  film.  Corning  sharp  cut  3-70  and  3-72 
filters  were  used  as  well  as  a  narrow -band  filter  («  50  A /half -width, 

X  egjr  **  5015  A).  Transmission  curves  for  the  3-70  and  3-72  filters 
are  shown  in  Fig.  21a  and  the  narrow -band  filter  curve  is  shown  in 
Fig.  21b. 


SECTION  IV 

DATA  ACQUISITION  AND  REDUCTION  TECHNIQUES 


All  test  parameters  were  scanned  in  approximately  0.  5  sec,  thus 
providing  essentially  simultaneous  reading  of  tunnel  and  plume  condi¬ 
tions.  These  data  were  digitized,  recorded  on  paper  tape,  and  subse¬ 
quently  reduced  on  the  VKF  CDC  1604 -B  computer. 


4.1  TUNNEL  CONDITIONS 

Test  section  flow  conditions  were  determined  from  free-stream 
impact  pressures  and  total  pressure  data  by  use  of  the  normal  shock 
relations.  Pressure  readings  were  correlated  for  orifice  and  probe 
effects  as  discussed  in  Appendix  I.  For  Txmnel  D,  total  temperature 
was  measured  with  a  thermocouple,  and  for  Tunnel  M  was  determined 
from  measurements  of  mass  flow  rate,  reservoir  pressure,  and  geo- 
metric  area  of  the  sonic  throat.  Using  calculations  similar  to  those 
reprorted  in  Ref.  8,  the  flow  was  assumed  to  be  in  equilibrium  upstream 
of  the  sonic  throat  and  frozen  downstream  with  a  specific  heat  ratio 
of  1.4.  -  -  - 
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Fig.  21  Transmission  Characteristics  of  Filters 
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4.2  PLUME  SOURCE  CONDITIONS 

Total  temperature  in  the  rocket  simulator  chamber  waa  measured 
by  means  of  a  thermocouple  which  passed  through  the  pressure  line  and 
extended  into  the  chamber  (Table  in).  Chamber  pressure  was  measured 
by  means  of  a  pressure  tap  in  the  chamber  wall.  For  the  0.  030-  and 
0.  080 -in. -throat -diam  simulators,  reservoir  pressure  could  be  taken 
equal  to  the  static  pressure.  However,  in  the  case  of  the  0.  3 15-in. - 
throat-diam  simtilator  with  a  chamber -to -thro  at  area  ratio  of  1.61,  the 
total-to- static  pressxore  ratio  was  calculated  to  be  1. 12  for  the  helium 
and  argon  involved  and  this  correction  was  applied  to  the  static  meas¬ 
urement. 

Average,  plus  maximTim  and  minimum,  values  of  total  tempera¬ 
ture  and  reservoir  pressure  in  the  simulator  chamber  are  given  in 
Table  IV  for  the  respective  conditions. 


4,3  PLUME  IMPACT  PRESSURE  SURVEYS 

An  axial  centerline  impact  pressure  survey  of  the  plume  was  made 
for  each  condition  with  measimements  made  at  l-in,  intervals  from  the 
simxilator  nozzle  exit  plane  (x  =  0)  to  the  Mach  disk.  Radial  surveys 
extending  from  the  centerline  (R  =  0)  past  the  outer  edge  of  the  test  core 
were  made  at  x  =  3  in  and  at  axial  stations  5  in.  apart  from  x  =  5  in.  to 
the  Mach  disk.  For  the  IC  conditions,  radial  surveys  were  also  made 
at  X  =  1  in.  In  the  mixing  region  of  the  plume,  measurements  were 
generally  made  at  nominal  1/8-in.  radial  intervals  with  1/4-  and  1/2-in. 
increments  being  used  in  regions  of  less  steep  gradients.  During  these 
surveys  probe  position  was  monitored  on  a  digital  voltmeter,  and  the 
probe  was  generally  within  ±0. 025  in.  of  the  nominal  position.  The 
nominal  radial  probe  position  was  used  in  the  tabulation  and  graphical 
presentation  of  impact  pressure  data  for  runs  1024  (condition  1C5)  and 
1025  (condition  IC2)  and  the  4000-  and  5000 -series  of  runs  (conditions 
IB  and  lA,  respectively)  because  excessive  noise  on  that  particular 
amplifier  of  the  analog-to-digital  converter  was  found  to  have  rendered 
the  recorded  signal  erratic. 

Since  the  properties  of  the  gas  mixture  in  the  plume  mixing  region 
were  unknown,  measmred  impact  pressures  were  corrected  for  viscous 
and  thermal  effects  as  explaned  in  Appendix  I,  assuming  first  that  only 
nitrogen  was  present  and  then  that  023ly  the  plume  gas  was  present. 
Although  the  former  could  be  correct  for  points  to  the  outside  of  the 
mixing  region  and  the  latter  for  points  to  fhe  inside,  within  the  mixing 
region  the  correct  value  should  fall  between  the  two  values  given. 


49 


TABLE  IV 
TEST  SUMMARY 


O 


I 


j 

I 


Chamber  Conditions 


Plume  Group  - 


Condition  Gas 

Series 

Test  Activity 

avg 

max 

min 

avg 

max 

min 

1C5 

Helium 

1000 

Impact  Pressure  Surveys 

764 

773 

753 

5.  04 

5.  30 

4.  93 

1500 

Flow  Visualization  . 

777 

791 

762 

4.  95 

5.  20 

4.  78 

1700 

Spectrometer  Measurements 

758 

770 

744 

4.  96 

5.  20 

4.  76 

LAB 

CO2 

5000 

Impact  Pressure  Surveys 

741 

750 

733 

396 

412 

380 

5500 

Flow  Visualization 

729 

740 

711 

402 

411 

386 

IIB 

Helium 

2 

Impact  Pressure  Sujrveys 

304 

305 

303 

55.  8 

56.  0 

55.  6 

IC2 

Helium 

1000 

Impact  Pressure  Surveys 

711 

722 

696 

2.  03 

2.  10 

1.91 

1500 

Flow  Visualization 

704 

721 

693 

1.  86 

2.  14 

1.  66 

1700 

Spectrometer  Measurements 

706 

723 

673 

1.  97 

2.  01 

1.  80 

lAC 

CO2 

5000 

Impact  Pressure  Surveys 

742 

756 

734 

97.  7 

107.  8 

89.2 

5500 

Flow  Visualization 

743 

749 

738 

97.  8 

99.4 

95.  5 

lAA 

CO2 

5000 

Impact  Pressure  Surveys 

731 

747 

716 

50.  5 

62.  9 

38.  0 

• 

5500 

Flow  Visualization 

731 

734 

728 

47.  6 

49.  1 

46.4 

IB 

Helium 

4000 

Impact  Pressure  Surveys 

723 

725 

720 

394 

394 

392 

4500 

Plow  Visualization 

719 

722 

717 

395 

396 

394 

4700 

Spectrometer  Measurements 

713 

721 

691 

398 

399 

396 

Argon 

3500 

Flow  Visualization 

551 

553 

547 

2,  94 

3.  16 

2.  80 

Flow  Visualization 

636 

641 

632 

7.  06 

7.  48 

6.  81 

Flow  Visualization 

725 

728 

721 

8.  00 

8.04 

7. 95 

1  i 


I 

I 

I 
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4.4  ELECTRON  BEAM 


When  making  flow  visualization  photographs  the  tunnel  and  plume 
generator  conditions  were  set  at  their  desired  values,  and  the  electron 
gun  housing  was  axially  traversed  until  the  electron  beam  was  located 
in  the  plume  generator  nozzle  exit  plane.  The  digital  voltmeter  indi¬ 
cating  the  axial  position  of  the  beam  was  set  to  zero.  Simultaneously^ 
the  camera  shutter  was  opened  and  the  electron  beam  was  started  on 
an  axial  traverse  downstream.  After  the  electron  beam  had  moved 
approximately  16  in.  downstream  from  the  plume  generator  nozzle 
exit,  the  axial  traverse  was  stopped  and  the  camera  shutter  was 
closed.  Tunnel  data  points  were  taken  throughout  the  scan  to  monitor 
tunnel  conditions  and  the  exact  length  of  axial  traverse. 

When  filters  were  used,  they  were  taped  directly  onto  the  camera 
lens  mount.  The  3-72  filter  was  normally  used  to  reduce  film  exposure 
by  N2"*’(1“)  and  N2(2'^)  predominantly  blue  radiation  yet  pass  the  pre¬ 
dominant  radiation  of  helium  at  5015-A  and  argon  at  4607-A  radiation. 
The  3-70  filter  was  used  to  block  all  radiation  below  4900  A.  The 
narrow-band  filter  passed  only  the  helium  5015-A  radiation;  when  it 
was  used,  several  beam  traverses  had  to  be  made  to  achieve  adequate 
film  exposure. 


For  rotational  temperature  measurements  the  spectrometer  slits 
were  set  at  20  microns  wide  by  10  mm  high.  The  rotational  spectra 
scans  were  made  at  a  scan  speed  of  2.  5  A/ min  with  a  lock-in  time  con¬ 
stant  of  1  sec  and  using  the  2-  to  20-mv  lock-in  scales.  For  a  given 
axial  position  and  plume  condition  no  more  than  two  rotational  spectra 
scans  were  made  at  a  radial  position. 


A  typical  rotational  spectra  is  shown  in  Fig.  22  (note  that  the 
R  branch  of  the  N2‘*’(l*')(0, 0)  band  was  always  used).  From  such  a 
scan  the  peak  heights  of  the  spectral  lines  are  measured.  The  peak 
heights  are  just  the  relative  intensity  values  of  the  rotational  lines 


(lV',K' 


,  );  using  these  values  and  the  known  values  of  , 

“A  V  ',  K  -1 


Eq.  (II-6)  of  Appendix  II  may  be  solved  for  rotational  temperature  by  an 
iterative  machine  calculation.  The  data  deck  requires  input  data  cards 
for  the  following:  spectral  scan  number  and  designation,  total  number 
of  spectral  lines  in  a  scan,  and  each  line's  relative  intensity  value. 
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The  intensities  of  even  numbered  rotational  lines  (K"  =  2,  4,  .  .  . )  are 
doubled  before  input  to  properly  take  account  of  the  influence  of 
nuclear  statistics.  Letting 


H  K' 


and  X  =  K  (K  +1),  the  machine  calculation  performs  a  least-squares 
fit  of  the  y  values  to  the  straight  line  given  by  Eq.  (II-6)(Fig.  23).  As 
shown  from  Eq.  (11-6)  the  slope  of  the  stright  line  is 

he  Bo 

"k  Tb 

i 

hence,  a  rotational  temperature  can  be  calculated  from  the  slope  of 
the  straight  line,  A  temperature  of  200®K  is  always  assumed  in  the 
first  step  of  the  iteration  process.  Using  200°K,  the  values  of  [G]q 
may  be  calculated  and  a  straight  line  is  then  fit  to  the  y  values.  From 
the  slope  of  the  line,  a  temperature  is  calculated.  If  this  calculated 
value  of  rotational  temperature  is  not  within  0.  01  percent  of  the 
assumed  value  of  rotational  temperature,  then  the  calculated  value  is 
used  as  the  next  assumed  value  of  temperature  and  the  calculation  is 
repeated. 

As  observed  by  Williams  (Ref.  II-9)  and  Ashkenas  (Ref.  9)  only 
for  temperatures  near  200*^  (pure  gas,  equilibrium  conditions)  is  the 
plot  of  y  values  versus  x  =K'  (K'  +1)  a  good  straight  linej  that  is,  the 
temperature  determined  from  the  slope  of  the  line  is  a  function  of  the 
number  of  rotational  lines  used.  For  this  reason  a  temperature  de¬ 
fined  by  the  maximum  niimber  of  spectral  lines  is  first  calculated  by 
the  iterative  procedure  described  in  the  preceding  paragraph.  Then 
the  entire  process  is  repeated  for  lines,  and  this  gives 

another  rotational  temperature.  This  process  is  repeated  i  times 
through  -i  +  1  =  6.  Only  if  the  variation  of  temperature  with  the 

number  of  rotational  lines  used  is  fairly  small  can  a  single,  well- 
defined  temperature  be  obtained. 

As  shown  in  Appendix  II  (Eq.  (II-8)),  a  calculation  of  band  intensity 
ratios  as  a  function  of  vibrational  temperature  can  be  made.  With 
these  calculations  a  measurement  of  band  intensity  ratios  will  yield 
the  vibrational  temperature.  For  band  intensity  measurements  the 
dove  prism  was  placed  in  the  optical  system,  and  the  spectrometer 
slits  were  normally  set  100  to  400  microns  wide  by  20  mm  high.  The 
output  of  the  lock-in  amplifier  was  fed  into  the  integrator  (Fig.  19), 
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and  the  output  of  the  integrator  was  displayed  by  a  digital  voltmeter. 
Before  a  scan  of  a  vibrational  band  the  light  from  the  electron  beam 
was  blocked  off,  and  the  integrator  null  was  adjusted  such  that  the 
integrator  output  was  constant  in  time.  The  integrator  output  was  then 
set  to  zero  with  the  "kill"  switch  (Fig.  19).  The  spectrometer  was 
started  in  its  scan  of  the  band  (at  a  rate  of  25  A/ min),  and  the  integrator 
was  placed  into  the  integrate  mode.  The  display  of  the  digital  voltmeter 
was  observed,  and  the  reading  was  recorded  after  the  band  was  com¬ 
pletely  scanned.  Beam  current  was  kept  constant  during  the  integra¬ 
tion  of  the  two  bands  to  be  used  for  the  intensity  ratio. 

The  nitrogen  number  density  was  determined  by  integration  of  the 
N2**"(1'')(0,  0)  band.  The  integration  technique  was  the  same  as  for  vibra¬ 
tional  temperature  determination.  Spectrometer  slits  were  60  microns 
wide  by  20  mm  high,  and  the  dove  prism  was  used.  At  the  end  of  each 
(0, 0)  band  integration  the  beam  current  was  measured  at  both  the  flat- 
plate  collector  on  the  timnel  floor  and  at  the  auxiliary  collector  (extended 
position)  located  directly  beneath  the  gun  orifice.  The  ratio  of  band  in¬ 
tensity  to  beam  current  (auxiliary  collector  value  used)  then  gave  a 
relative  value  of  nitrogen  number  density,  as  shown  by  Eq.  (11-13), 
assuming  Ty  to  have  the  same  value  in  the  plume  region  as  in  the  free- 
stream  region.  Absolute  nitrogen  number  densities  are  then  obtained 
by  relating  these  relative  values  to  the  band  intensity-beam  current 
ratio  obtained  from  centerline  free-stream  measurements  and  the 
Timnel  M  data  reduction  program  calculation  of  free-stream  nitrogen 
density  (Appendix  II). 

The  determination  of  helium  number  density  was  obtained  by  the 
integration  of  the  He  50 15 -A  line.  The  integration  technique  was  the 
same  as  that  for  nitrogen  number  density  determination.  Spectrometer 
slits  were  60  microns  wide  by  20  mm  high,  and  the  dove  prism  was 
used.  Again,  beam  current  was  measured  at  both  positions  at  the  end 
of  each  integration.  The  ratio  of  line  intensity  to  beam  current 
(auxiliary  collector  value  used)  then  gave  a  relative  value  of  heliinn 
number  density,  as  shown  by  Eq,  (II-l).  Absolute  helium  number  den¬ 
sity  values  were  then  obtained  using  the  calibration  of  the  ratio  of  line 
intensity  to  beam  current  versus  helium  number  density  (Appendix  ID. 

SECTION  V 

EXPERIMENTAL  RESULTS 


5,1  PHOTOGRAPHS 

Color  still  photographs  were  taken  for  all  conditions  in  Tunnel  M 
using  the  electron  beam  as  the  source  of  excitation.  PTiotographs  were 
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taken  with  and  without  filters,  as  discussed  in  Section  4,4,  Photo¬ 
graphs  taken  with  the  filters  showed  considerably  more  detail  than 
those  without  filters.  Photographs  taken  with  the  various  filters  were 
usually  quite  similar,  differing  only  in  details  of  coloring  which  could 
be  altered  in  processing.  As  far  as  revealing  the  basic  structure  of 
the  flow  field  in  the  photography,  the  specific  filter  used  was  found  to 
be  of  minor  importance. 

Selected  photographs  are  given  in  Figs.  24  through  28.  The  first 
three  of  these  are  for  helium  plumes  and  are  for  conditions  ICS,  IC2, 
and  IB.  The  fourth  is  for  a  CO2  plume,  condition  lAB,  and  the  fifth  is 
of  an  argon  plume. 

Figure  24  gives  the  photograph  for  Group  1517,  condition  ICS, 
using  the  3-70  filter.  In  the  coloring  of  this  particular  photograph 
nitrogen  appears  as  red  and  helium  as  green,.  The  outer  shock 
appears  quite  prominently  as  the  line  dividing  the  bright  red  region 
from  the  background.  The  nitrogen  shock  layer  appears  as  the  red 
region.  In  the  region  near  the  nozzle  exhaust,  the  inner  shock  appears 
as  the  line  dividing  the  bright  green  region  from  the  imer  dark  green 
background.  The  dividing  line  between  the  nitrogen  and  helium  is  quite 
distinct. 

As  one  moves  farther  downstream,  the  inner  shock  becomes  less 
distinct  because  the  density  is  low.  It  is  not  clear  whether  the  shock  is 
becoming  thicker  as  one  moves  downstream  or  the  densities  are  so  low 
that  the  photographs  do  not  register  the  small  difference  in  density 
across  the  shock.  The  contact  region,  which  begins  as  a  thin  line  near 
the  exit  plane,  becomes  very  thick.  The  outer  limits  of  the  helivun 
and  the  inner  limit  of  the  nitrogen  can  be  picked  off  the  photograph  as 
the  limit  of  the  appropriate  color. 

The  large  thickness  of  the  mixing  region  would  not  be  anticipated 
on  the  basis  of  experience  under  lower  altitude  conditions.  Farther 
downstream,  the  thickness  of  the  mixing  region  is  several  times  that 
of  the  nitrogen  shock  layer. 

Figure  25  is  a  photograph  for  condition  IC2,  which  is  identical  to 
condition  ICS  except  for  the  chamber  pressure,  which  is  2  psi  rather 
than  5  psi.  Although  the  photograph  is  in  every  way  similar,  a  large 
relative  plume-length-is  shown.- -The  closing  of  the“inner“Shock  is  faintly 
visible. 

Figiire  26  is  for  condition  IB,  which  is  for  a  high  chamber  pressure 
(450  psi)  and  the  large  area  ratio  nozzle.  The  inner  shock  is  consider¬ 
ably  more  distinct  than  in  the  two  previous  photographs  because  of  the 
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16  in. 


Fig.  26  Photograph  of  Plume  for  Condition  IB.  Tunnel  M.  He  Plume, 
A,/A*  =  25.  d*  =  0.030  in.,  and  p..  *  450  psia 


Fig.  27  Photograph  of  Plume  for  Condition  lAB,  Tunnel  M.  CO2  Plume, 
A,/A*  =  25.  d*  =  0.030  in.,  and  p^  =  400  psia 
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Fig.  28  Photograph  of  Argon  Plume,  Tunnel  M,  A^/A 
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considerably  higher  density.  The  mixing  region  is  stHl  relatively 
large.  The  long  slender  nature  of  this  plume  is  appropriate  for  the 
large  value  of  X  for  this  particular  combination  of  nozzle  area  ratio 
and  7. 

Figure  27  is  for  the  condition  lAB,  which  is  a  CO2  plume  with 
conditions  chosen  to  match  condition  ICS.  Because  of  the  differences 
in  d  and  p^,,  direct  comparison  of  physical  lengths  on  the  correspond¬ 
ing  photographs  (Figs.  24  through  27)  should  not  be  made. 

Figure  28  is  for  a  condition  similar  to  IC2  except  with  argon  as 
the  plume  gas.  This  particular  condition  was  run  in  an  unsuccessful 
attempt  to  evaluate  effects  of  condensation.  The  photograph  is  included 
because  it  shows  a  larger  extent  of  the  inner  plume  than  any  of  the 
others.  In  particular,  regions  around  the  Mach  disk  are  shown.  It  is 
not  possible,  however,  to  pick  out  a  distinct  line  corresponding  to  the 
Mach  disk, 

No  worthwhile  photographic  data  were  retained  from  Tunnel  D. 
Color  schlieren  photographs  were  taken,  but  were  too  faint  to  be  of 
value.  The  installation  photograph  (Fig.  13)  was  actually  taken  during 
a  run,  and  it  is  obvious  that  little  flow  information  can  he  drawn  from 
the  photograph. 


5.2  IMPACT  PRESSURE  SURVEYS 

The  corrected  impact  pressure  data  for  the  various  conditions  are 
given  in  the  series  of  figures  from  Fig.  29  to  Fig.  35.  For  the  data 
taken  in  Tunnel  M,  the  pressure  data  have  been  corrected  for  orifice 
and  probe  effects  (Appendix  I).  In  the  outer  flow  the  raw  data  have 
been  corrected  using  values  for  nitrogen,  whereas,  for  the  inner  flow, 
values  appropriate  for  the  plume  gas  have  been  used.  In  the  mixing 
region  values  are  given  for  both  gases. 

For  each  condition  tested  in  Ttmnel  M  the  (a)  part  of  the  figure 
gives  an  axial  centerline  traverse,  and  the  succeeding  parts  give  results 
for  radial  surveys  at  various  axial  locations. 

For  the  axial  surveys  the  most  striking  feature  is  the  minimum  in 
the  value  of  impact  pressure.  This  minimum  is  assumed  to  occur  at 
the  location  of  the  Mach  disk.  Table  V  summarizes  the  results  for  the 
location  of  the  Mach  disk.  Also  shown  is  the  Mach  disk  location  normal¬ 
ized  with  the  scaling  parameters  of  Eq.  (19).  For  the  first  condition 
of  Table  I  (conditions  ICS  and  lAB),  the  agreement  is  adequate,  but  for 
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a.  Axial  Centeriine  Corrected  Impact  Pressure  Survey 
Fig.  29  Impact  Pressure  Surveys  for  Condition  ICS 
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p^,  mm  Hg 

b.  Radial  Corrected  Impact  Pressure  Survey  {x  =  1  in.) 


c.  Radial  Corrected  Impact  Pressure  Survey  (x  =  3  in.) 
Fig.  29  Continued 
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y.  in. 


Ouier 

Shock 


Inner 

Shock 


d.  Radial  Corrected  Impact  Pressure  Survey  (x  =  5  in.) 


p^',  ram  Hg 


e.  Radial  Corrected  Impact  Pressure -Survey  (x  =  10  in.)  - 
Fig.  29  Continued 
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f.  Radial  Corrected  Impact  Pressure  Survey  (x  =  15  in.) 
Fig.  29  Concluded 
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Fig.  30  Impact  Pressure  Surveys  for  Condition  lAB 
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b.  Radial  Corrected  Impact  Pressure  Survey  (x  =  3  in.) 
Fig.  30  Continued 
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p^.  mmHg 

d.  Radial  Corrected  Impact  Pressure  Survey  (x  °  10  in.} 
Fig.  30  Continued 
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a.  Radial  Impact  Pressure  Survey  <x  =  3  in.) 
Fig.  31  Impact  Pressure  Surveys  for  Condition  MB 
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b-  Radial  Impact  Pressure  Survey  (x  =  5  in.) 
Fig.  31  Continued 
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c.  Radial  Impact  Pi«ssure  Survey  (x  -  10  in.) 
Fig.  31  Conduded 
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a.  Axial  Centerline  Corrected  Impact  Pressure  Survey 
Fig.  32  Impact  Pressure  Surveys  for  Condition  IC2 
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c.  Radial  Corrected  Impact  Pressure  Survey  (x  =  3  in.) 
Fig.  32  Continued 
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a.  Axial  Centsriine  Corrected  Impact  Pressure  Survey 
Fig.  33  Impact  Pressure  Surveys  for  Condition  lAC 


76 


AeDC-TR-71-25 


b.  Radial  Corrected  Impact  Pressure  Survey  (x  =  3  in.) 
Fig,  33  Continued 
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c.  Radial  Corrected  Impact  Pressure  Survey  (x  =  5  in.) 
Fig.  33  Continued 
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e.  Radial  Corrected  impact  Pressure  Survey  (x  ==  15  in.) 
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8.  Axial  Centerline  Corrected  Impact  Pressure  Survey 
Fig.  34  Impact  Pressure  Surveys  for  Condition  lAA 


8 


AEDC-TR-71-25 


b.  Racial  Correctad  Impact  Pressure  Survey  (x  =  3  in.) 
Fig.  34  Continued 
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c.  Radial  Corrected  Impact  Pressure  Survey  (x  =  5  in.) 
Fig.  34  Continued 
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d.  Radial  Corrected  Impact  FVessure  Survey  (x  =  10  in.) 
Fig.  34  Continued 
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e.  Radial  Corrected  Impact  Pressure  Survey  {x  =  12  in.) 
Fig.  34  Concluded 


AEDC-TB-71-25 


c.  Ratfial  Corrected  impact  Pressure  Survey  (x  =  5  in.) 
Rg.  35  Continued 
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d.  Radial  Corrected  Impact  Pressure  Survey  (x  10  in.) 


e.  Radial  Corrected  Impact  Pressure  Survey  (x  =  15  in.) 
Rg.  35  Concluded 
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TABLE  V 

EXPERIMENTAL  DATA  ON  MACH  DISK  LOCATION 
AS  DETERMINED  BY  IMPACT  PRESSURE  SURVEYS 


Condition 

d*,  in. 

X  Mach 
Disk, 
in. 

A 

X 

r*  % 

1  Pc  ^Fm 

ICS 

0.  315 

15.0 

4.  5 

5.  55 

lAB 

0.030 

14.0 

4.  7 

5.  41 

IC2 

0.  315 

8.  0 

4.  5 

4.  70 

lAC 

0.  080 

17.  0 

4.  7 

4.  95 

lAA 

0.  080 

11.  5 

4.7 

4.  74 

IB 

0.  030 

18.  0 

33.  5 

7.  79 

S9 
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the  second  base  condition  of  Table  I  (conditions  IC2,  lAC,  and  lAA), 
the  agreement  is  reasonable.  Although  there  are  differences  in  the 
values  for  the  two  sets  of  data  for  X***  4.  6,  the  average  of  all  values 
of  (x/ r*)  y}{qj (Pc  at  4.  6  is  5. 07  ±  0.42,  About  half  of  this 

spread  could  be  caused  by  variations  in  effective  r*  attributable  to  low 
Reynolds  numbers  in  the  rocket  simulator  throats,  and  the  remainder 
could  be  brought  about  by  other  nonuniform  flow  conditions  in  the  small 
nozzles.  External  low  density  flow  phenomena  caimot  be  ruled  out 
either. 

The  radial  pressure  surveys  of  Fig.  29  are  typical  of  the  other 
surveys  in  terms  of  gross  details.  Figure  29b  gives  data  taken  in  the 
plane  1  in.  from  the  exit  plane.  In  this  figure  the  prominent  features 
are  the  inner  shock,  the  outer  shock,  and  the  timnel  shock.  These  can 
be  noted  in  the  photograph  of  Pig.  24.  The  (c)  part  of  the  figtire  gives 
simitar  results,  with  the  inner  shock  and  the  outer  shock  displaced 
farther  from  centerline  and  the  tunnel  shock  closer.  There  is  an  indi¬ 
cation  of  what  could  be  the  contact  surface.  In  the  (d)  part  of  the  fig¬ 
ure,  there  is  no  indication  of  the  contact  sttrface.  In  Fig.  29e  the 
inner  shock  appears  as  a  smeared  region,  rather  than  a  discontinuity. 
The  tunnel  shock  and  the  outer  shock  have  merged.  This  merger  does 
not  appear  at  this  location  in  the  photograph  (Fig.  24)  because  the 
tunnel  shock  intersects  the  centerline  nearer  the  exit  with  the  probe  in 
the  tunnel  because  of  a  higher  tank  presstire  caused  by  the  probe  block¬ 
age  effect.  In  the  (f)  part  of  the  figure  everything  has  smeared  out. 

This  location  is  dose  to  the  Mach  disk  location  (Fig.  29a).  The  results 
for  the  remaining  figures  are,  in  general,  similar  to  those  of  Fig.  29. 


5.3  DENSITY  AND  TEMPERATURE  MEASUREMENTS 
5.3.1  Density 

Using  the  electron  beam  flourescence  technique  in  Tunnel  M,  the 
radial  profiles  of  the  number  density  of  neutral  ground  state  atomic 
helium  and  molecular  nitrogen  were  determined  at  the  axial  positions 
located  3,  5,  and  10  in.  downstream  of  the  exit  plane  of  the  plume 
nozzle  for  plume  conditions  IC2  and  ICS,  and  similar  data  were  obtained 
3  in.  from  the  exit  plane  of  the  plume  nozzle  for  condition  IB, 

Absolute  number  density  data  were  obtained  in  the  manner  described 
in  Appendix  IE  with  the  experimental  apparatus  discussed  in  Section  3.  5, 
using  the  raw  intensity  data  of  the  atomic  and  molecular  radiation  and 
measured  beam  current.  Figures  36a,  b,  and  c  present  the  radial  pro-  ■ 
files  of  the  absolxrte  ntimber  densities  for  condition  ICS;  Figs.  37a,  b, 
and  c  for  condition  IC2;  and  Fig.  38  for  condition  IB.  The  data  given  in 
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O  N2  Number  Density 
#  He  Number  Density 


Condition  1C5 
Tunnel  M 
He  Plume 
Ae/A*  •  1.61 
d*  •  0.  315  In. 
•  5  psia 


a.  X  =  3  in. 

Fig.  36  Electron  Beam  Number  Densities  for  Condition  ICS 
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b.  X  =  5  in. 
Fig.  36  Continued 
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•  5  pei* 
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c.  X  =  10  in. 
Fig.  36  Concluded 


a.  X  =  3  in. 

Fig.  37  Electron  Beam  Number  Densities  for  Condition  IC2 
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c.  X  =  10  in. 
Rg.  37  Concluded 
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Fig.  38  Electron 
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these  figures  were  obtained  assuming  the  line  or  band  intensity  to  be  a 
linear  function  of  number  density;  i.  e. ,  no  allowance  has  been  made 
for  coUisional  quenching  effects.  Figures  39  through  41  exhibit  the 
same  experimental  data  where  coUisional  quenching  effects  are  esti¬ 
mated  using  the  quenching  ratio  constants  listed  in  Table  Il-in, 

Appendix  II.  It  should  be  noted  that  until  low  temperature  rate  con¬ 
stants  are  e:q)erimentally  determined.  Figs.  39  through  41  indicate 
only  an  estimate  of  the  effect  of  collisions. 

Listed  in  Table  VI  and  shown  in  Fig.  42  are  the  locations  of  the 
maxima  of  the  number  density  in  the  radial  direction  as  a  function  of 
axial  position  for  all  plume  conditions  investigated. 

Figures  43a  and  b  present  the  centerline  values  of  helium  number 
density  as  a  fxmction  of  the  inverse  square  of  the  axial  position. 

5.3.2  Rotational  Temperature 

Radial  profiles  of  the  rotational  temperatures  of  the  ground  state 
neutral  nitrogen  molecule  were  determined  in  Tunnel  M  using  the  (0, 0) 
band  of  the  first  negative  system  of  N2'**  as  described  in  Appendix  II  for 
the  IC2  and  ICS  conditions  at  the  axial  position  10  in.  downstream  from 
the  plume  nozzle  exit  plane  and  for  the  ICS  condition  at  the  5-in.  axial 
position.  These  data  are  shown  in  Figs.  44a,  b,  and  c,  where  the 
rotational  temperatures  obtained  from  iterated  computer  calculations 
described  in  Appendix  n  are  plotted  as  a  function  of  the  number  of 
spectral  lines  employed  for  the  determination.  If  a  rotational  tempera¬ 
ture,  Tjj,  is  defined,  the  Tp  variation  with  K'  should  exhibit  a  zero 
slope^  whereas  a  finite  slope  is  indicative  of  the  non-Boltzmann  nature 
of  the  radiating  species  and  can,  in  all  probability,  be  attributed  to  both 
secondary  electrons  and  coUisional  processes  of  extremely  large  cross 
sections.  Since  these  cross  sections,  or  rate  constants,  are  presently 
unknown,  no  density  corrections  have  been  applied  to  the  data. 

The  rotational  temperatiire  obtained  using  ten  spectral  lines,  serv¬ 
ing  only  to  show  the  gross  character  of  rotational  energy  distribution 
throughout  the  plume,  is  shown  in  Fig.  45  as  a  fimction  of  the  radial 
position  for  the  IC2  and  ICS  conditions  at  the  10 -in.  axial  position,  and 
these  values,  as  weU  as  the  ICS  data  at  5  in. ,  are  listed  in  Table  VII. 

5.3.3  Vibrational  Temperature 

Using  the  (0, 1)  and  (1, 2)  bands  of  the  first  negative  system  of  N2**’, 
electronically  integrated  vibrational  band  intensity  ratios  were  obtamed 
for  both  the  IC2  and  ICS  conditions  at  the  3-  and  5 -in.  axial  positions. 
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a.  X  =  3  in. 

Fig.  39  Eiectron  Beam  Number  Densities  with  Coilisional  Quenching  Effects  for  Concfition  tC5 
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b.  X  =•  5  In. 
Fig.  39  Continued 
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c.  X 
Fig.  39 
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a.  X  3  in. 

Rg.  40  Electron  Beam  Number  Densities  with  Collisional  Quenching  Effects  for  Condition  1C2 
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Fig.  40 
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c.  X  =  10  in. 
Fig.  40  Concluded 
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O  N2  Number  Density,  Corrected 
#  He  Number  Density,  Corrected 
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Tunnel  M 
He  Plume 
A^/A*  =  25.  0 

d*  =  0.  080  in. 
=  450  psia 


Fig.  41  Electron  Beam  Number  Densities  with  Collisional  Quenching  Effects  for  Condition  IB  (x  =  3  in.l 
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Fig.  42  Locations  of  tfie  Maxima  of  the  He  and  N2  Number  Densities 
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TABLE  V! 

ELECTRON  BEAM  RESULTS 


Location  of  Peaks  y,  in. 


IC2 

IC5 

IB 

X,  in. 

N2  He 

Ng  He 

Ng  He 

3 

1.8  1.2 

1.8  1.0 

1.6  0,9 

5 

2.4  1.4 

2.  6  1.  8 

10 

3.5  1.6 

4.1  2.2 

107- 


N  X  10”^^,  molecules/c 
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a.  Condition  ICS 

Fig.  43  Centerline  Number  Density  of  He  versus  Inverse  Square  of  Axial  Distance 
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b.  Condition  IC2 
Fig.  43  Concluded 
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Fig.  45  N2  Rotational  Temperature  versus  Radial  Position  Using 
Ten  Spectral  Lines 
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TABLE  VII 

TEN  SPECTRAL  LINE  Tr  AS  FUNCTION  OF  RADIAL  POSITION 


Plume  Condition: 

ICS 

Axial  Position: 

5  in. 

y.  in. 

1.  97 

Tr,  ok 

185 

Plume  Condition: 

IC2 

Axial  Position: 

10  in. 

y.  in. 

2.  37 

2.76 

3.  16 

3.55 

3.  95  4.  35 

Tr,  ok 

200 

214 

172 

102 

85  80 

Plume  Conditions; 

ICS 

Axial  Position: 

10  in. 

y,  in. 

3. 16 

3.55 

3.  95 

4.  74 

Tr,  ok 

209 

250 

177 

99 

114 
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The  band  intensity  ratio  yielded  the  vibrational  temperature  by  use  of 
Fig.  II-6  of  Appendix  11.  Density  corrections  based  on  rate  constants 
experimentally  obtained  at  300“K  and  estimated  when  unknown  show 
that  the  density  correction  to  is  approximately  3  percent;  the 
smallness  of  the  correction  is  caused  primarily  by  the  assumed 
equality  of  ki(N2'^  j  He)  and  kg (N2'*',  He),  defined  in  Appendix  H,  and 
this  assumption  is  made  solely  because  of  ignorance  of  the  correct 
values.  Consequently,  no  density  correction  has  been  applied  to  the 
data. 

Table  VIII  lists  the  average  band  intensity  ratios  and  the  corre¬ 
sponding  vibrational  temperatures  as  a  function  of  the  radial  position 
for  the  various  conditions  investigated,  and  these  residts  are  presented 
injPigs,  46  and  47.  Also  listed  in  Table  VII  is  the  parameter  AT^ 
which  is  a  measure  of  the  re peatibility  of  the  result  when  more  than  one 
value  of  the  ratio  was  obtained. 


5.4  COMPARISON  OF  RESULTS 

Two  types  of  comparisor^  of  results  are  made.  The  first  is  the 
comparison  of  results  from  different  techniques  for  the  same  condition, 
and  the  second  is  comparison  of  results  for  different  conditions  in  terms 
of  the  scaling  parameters. 

5A1  Results  frmti  Different  Techniques 

These  comparisons  are  made  of  the  tests  in  Tunnel  M  using  re¬ 
sults  from  impact  pressure  surveys,  density  measurements  by  the 
electron  beam  technique,  and  photographs.  Figure  48  gives  compari¬ 
sons  for  condition  IC,  Figure  48a  is  at  the  x  location  of  3  in.  Impact 
pressures  are  from  Fig.  29c,  and  the  helium  and  nitrogen  densities 
are  from  Fig.  36a.  The  photographic  results  from  Fig.  24  are  shown 
as  bars  -  a  solid  bar  for  bright  and  a  dashed  one  for  what  appears  as 
the  mixing  region.  Figure  48b  is  f or  x  =  5  in.  and  Fig.  48c  for  10  in. 

It  appears  as  if  the  peak  in  the  impact  pressure,  the  peak  in  the 
nitrogen  density,  and  the  outer  edge  of  the  bright  brown  region  are  at 
about  the  same  distance  from  the  centerline.  This  location  may  then 
be  properly  identified  with  the  outer  edge  of  the  shock. 

The  peak  in  the  heliimi  density  does  not  seem  to  coincide  exactly 
with  any  other  measurement.  Figure  48d  summarizes  the  results  for 
condition  IC5,  Figs.  49a  through  d  are  analogous  to  Figs.  ■48a' through  d 
for  IC2,  and  Fig.  50  is  for  condition  IB  at  x  =  3  cm.  Photographic 
information  has  not  been  included  on  this  figure  because  the  plume  is  so 
thin  (Fig.  28). 
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Rg.  46  N2  Vibrational  Temperature  versus  Radial  Position  for 
Condtion  IC2 


Fig.  47  N2  Vibrational  Temperature  versus  Radial  Position  for 
Condition  ICS 
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TABLE  VIII 

VIBRATIONAL  TEMPERATURE  AS  A  FUNCTION  OF  RADIAL  POSITION 

Plume  Condition  IC2 


Axial  Position,  in. 


Radial 

3 

5 

Position 
y,  in. 

Ratio 

Ty.  “K 

AT^,  “K 

Ratio 

Ty,  “K 

ATv. 

“K 

1.  18 

2,  06 

3320 

60 

— 

— 

— 

1.  58 

2.  37 

2780 

40 

2.  08 

3440 

— 

1.97 

2.42 

2700 

20 

2.37 

2780 

— 

2.37 

2.  32 

2860 

40 

2.37 

2780 

— 

2.76 

— 

— 

— 

2.42 

2700 

40 

Plume  Condition  IC5 

Axial  Position,  in. 

Radial 

3 

5 

Position 
y.  in. 

Ratio 

T^.  “K 

AT^.  «K 

Ratio 

Ty.  "K 

ATy, 

“K 

1.58 

2.  15 

3140 

0 

2.24 

3000 

60 

1.  97 

2.  37 

2780 

— 

— 

— 

— 

2.  37 

2.  38 

2760 

20 

2.  22 

3040 

20 

2.  77 

— 

— 

— 

2.  28 

2900 

40 

3.  16 

— 

— 

2.  20 

3050 

20 
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Fig.  48  Comparison  of  Results  for  Condition  ICS 
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Figure  51  gives  a  comparison  of  the  shock  locations  for  various 
conditions  and  gases,  with  the  distances  normalized  by  the  scale  length, 
L  (Eq,  (17)).  The  values  for  helium  seem  to  correlate  fairly  well,  as 
do  those  for  CO2.  However,  the  locations  for  the  CO2  are  about  20  per¬ 
cent  less  than  those  for  helium.  The  reason  for  this  has  not  been  ex¬ 
plained. 


x/n 

Fig.  51  Normalized  Location  of  Maximum  Impact  Pressure 


SECTION  VI 

SUMMARY  AND  CONCLUSIONS 


Scaling  parameters  for  rocket  exhaust  plumes  were  determined  on 
the  basis  of  a  simplified  inviscid  theory.  These  parameters  were  used 
to  determine  test  conditions  for  the  simulation  of  the  flight  of ’a  typical 
rocket  through  the  atmosphere  at  an  altitude  of  80  km  at  Mach  numbers 
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of  18  and  10.  For  the  test,  CO2  and  He  were  chosen  for  the  e:diaust 
gases  of  the  models;  CO2  was  chosen  because  its  properties  closely 
matched  those  of  the  actual  rocket  plume,  and  He  was  chosen  because 
of  the  suitability  for  density  and  temperature  surveys  with  the  eiect'ron 
beam. 

Models  with  three  different  geometries  and  with  the  two  exhaust 
gases  were  tested  in  two  different  wind  tunnels  (M  and  D).  Rocket 
chamber  conditions  also  were  varied  to  enable  determining  plume 
properties  under  a  variety  of  conditions.  Pitot  pressure  surveys  were 
made  at  all  test  conditions;  photographs  were  taken  in  Tunnel  M  using 
the  electron  beam  for  flow  visualization;  and  densities  and  tempera¬ 
tures  were  determined  during  the  helium  rocket  exhaust  test  in 
Tunnel  M  using  radiation  measurements  from  the  electron  beam 
excited  flow. 

The  shocks  were  found  to  be  somewhat  diffuse  and  thick,  and  it  is 
believed  that  there  were  rarefaction  effects  on  plume  size  and  shape 
(characteristics  solutions  for  the  exact  conditions  would  be  necessary 
to  establish  this).  Density  measurements  of  the  free-stream  and 
plume  gases  demonstrated  that  significant  mixing  of  the  two  gases  had 
occurred  across  the  interface.  The  test  produced  useful  detailed  data 
on  simulated  rocket  exhaust  plvimes  interacting  with  a  rarefied  hyper¬ 
sonic  gas  stream. 
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Under  rarefied  flow  conditions  impact  pressure  probes  may  not 
indicate  the  ideal  impact  pressure  because  of  viscous  and  transitional 
flow  phenomena.  The  viscous  effects  are  attributable  to  the  thickening 
of  the  boundary  layer  and  shock  and  their  eventual  merging.  At  lower 
Reynolds  nmnbers  the  shock  becomes  more  diffuse  and  the  flow  departs 
further  from  equilibrium  as  free  molecular  conditions  are  approached. 
In  the  present  report,  all  free  molecular  conditions  are  based  on  equa¬ 
tions  from  Ref.  I-l.  The  specific  effects  considered  are  the  external 
or  probe  effect  which  was  reported  by  Potter  and  Bailey  (Ref.  1-2),  and 
the  orifice  effect  discussed  first  by  Potter,  Kinsiow,  and  Boylan 
(Ref.  1-3)  for  a  nonequilibrium  state  attributable  to  a  flux  of  energy  at 
a  surface.  The  ideal  impact  pressure  is  needed  in  order  to  determine 
conditions  in  the  free  stream,  so  that  the  measured  values  must  be 
corrected  for  the  probe  and  orifice  effects. 

Various  papers  have  been  written  on  the  probe  effect,  showing 
first,  as  one  proceeds  to-^ard  more  rarefied  flow,  a  small  decrease 
in  the  measured  impact  pressure  below  the  ideal  and  then  a  larger 
pressure  rise  as  near-free  molecule  conditions  are  encountered. 

The  interpretation  of  impact  probe  measurements  is  further  com¬ 
plicated  by  the  orifice  effect.  This  effect  is  manifested  by  an  inequality 
between  the  pressure  sensed  by  an  orifice  and  the  pressure  existing  at 
the  surface  containing  the  orifice  in  the  presence  of  a  flux  of  energy  or 
momentum.  The  orifice  effect  was  reported  in  Ref.  1-3  for  the  case  of 
the  energy  flux  to  a  surface. 

The  method  used  to  correct  pitot  probe  measurements  for  the 
probe  and  the  orifice  effect  to  obtain  the  true  impact  pressure  follows. 

Probe  Effect 


A  study  of  impact  probes  under  conditions  similar  to  those  in 
Tunnel  M  was  reported  by  Daun,  et  al.  (Ref.  1-4).  An  empirical  fit 
to  the  Daun  data  and  to  data  from  Tunnel  M  resulted  in  the  following 
equations  for  flat-faced  probes. 


91.4 


pw 

A 


L —  +  IIP*  (pz-'Poa) 
(Pw/Po  >Im  ■  " 


■  (I-l) 


91.4 


',2  "  1-5  R':2,r  VPj/Pm  +  »<'2.rp 

{pw/Po  j  £m 
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in  which  is  the  pressure  at  the  face  of  the  probe,  obtained  from  the 
indicated  pressure  by  applying  the  orifice  correction,  and  Pq'is  the  ideal 
impact  pressure. 


The  free  molecule  limit  of  p^/pQ'for  large  Mach  numbers  is 


y±i.  _  - - 

/Pw\  /  2  \y-i  y-i/  /  B-  T,  1  \ 

(po/ft.  "(y  +  0  ^  \  \  2y  ^ 


Y_ 

y-] 


<I-2) 


In  order  to  evaluate  these  equations  an  iteration  scheme  was 
necessary  since  M„,  Re2,  and  P2lPm  require  the  knowledge  of 
Pq'.  The  numerical  solution  is  discussed  below. 


Orifice  Effect 


This  effect,  attributable  to  the  nonequilibrium  state  of  the  gas  at 
the  solid  boundary,  when  a  flux  of  energy  is  present,  causes  the  pres¬ 
sure  sensed  by  an  orifice  to  deviate  from  the  force  per  unit  area  on 
the  surface.  The  amount  of  deviation  is  expressed  in  terms  of  a 
dimensionless  pressure  parameter,  p,  defined  by 

P  =  -  P5£„)/(Pw  -  Pii^) 


where 

Pi  is  the  pressure  inside  the  orifice  cavity, 

^ifm  value  of  p^  for  a  free  molecule  orifice,  and 

is  the  true  wall  pressure  in  the  absence  of  the  orifice. 

In  Refs.  1-3  and  1-5,  p  was  shown  to  be  a  function  of  the  Knudsen  num¬ 
ber  based  on  the  orifice  diameter,  d,  the  wall  temperature,  T^,  and 
Pi*  However,  recent  measurements  of  the  orifice  effect  in  Tunnel  M 
and  data  from  Cornell  Aeronautical  Laboratory  (Ref.  1-6)  show  that  the 
results  of  the  correlation  parameter  used  in  Refs.  1-3  and  1-5  are 
inadequate  for  higher  enthalpy  flows.  The  method  of  correcting  for 
the  orifice  effect  in  this  test  is  based  upon  a  paper  presented  at  the 
Seventh  International  Symposiiim  on  Rarefied  Gas  Dynamics  at  Pisa, 
Italy,  June  30,  1970  (Ref.  1-7). 

The  relation  between  p  and  a  Knudsen  number,  Kn,  based  upon 
orifice  diameter,  p^,  and  a  reference  temperature  Tj,  that  was  used  is, 

P  =  Z/(l  +  Z)  (1-4) 
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in  which 


7  0.394  0.282  0.038 

z  =  — — —  +  - + 


Ks 


Kn^ 


and  Tj,  are  both  functions  of  k,  the  dimensionless  heat-transfer 

rate  parameter,  and  a,  the  thermal  accommodation  coefficient. 

Pifjjj/  Pw  and  Tp/T-w  are  given  by: 


Pw 


—  a) 


+  4 


(1-5) 


2  -  TTK^ 


and 


=  1  +  2j7/c^  ~  ^  “^2^  l2ffK*(l  -  a)  +  or] 


(1-6) 


where 


q  <y  ~  n 
P^{y  +  l)v  RTw 


(1-7) 


An  iteration  scheme  is  required  to  solve  for  p^.  The  thermal  accom¬ 
modation  coefficient,  a,  was  presumed  to  be  0.  79  for  N2,  0.  8  for  CO2 
and  0.  51  for  He.  The  values  for  helium  and  nitrogen  are  those  used 
in  Ref.  1-5  and  close  to  those  obtained  in  Ref.  1-8.  The  carbon  dioxide 
value  was  obtained  from  a  consideration  of  molecular  weight. 


Based  on  limited  heat-transfer  data  from  Tunnel  M  for  flat -faced 
bodies,  the  following  equation  was  used  to  calculate  the  stagnation 
point  heat  transfer 

qjm 

Hk^FR/ 

where  is  the  free  molecule  heat-transfer  rate  and  q  is  from  the 

continuum  Fay-Riddell  equation  (Ref.  1-9).  This  experession  gives  the 
correct  free  molecule  heat -transfer  rate  and  approaches  the  Fay-Riddell 
theory  for  higher  Reynolds  numbers. 


The  equations  used  for  q^^  and  q  pjj  for  higher  Mach  numbers  are 

*  7^  -  ^  VSJT  M„]  -  (1-9) 


qfm  =  -«Po 
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and  based  upon  Ref.  1-9,  the  heat -transfer  rate  to  a  flat-faced  probe  of 
radius  r  is  _ 


(I- 10) 


based  upon  the  results  of  Ref.  1-9.  Heat-transfer  rate  is  negative  for 
"cold  wall"  conditions. 

If  free -stream  properties  are  known  then  Eqs.  (I-l)  through  (I- 10) 
can  be  solved  to  permit  the  correction  of  data  for  the  orifice  and  probe 
effects.  However,  since  free-stream  conditions  are  known  only  after 
applying  the  correction,  an  iteration  scheme  was  used.  The  numerical 
method  for  correcting  the  pitot  data  is  described  briefly  below. 

Iteration  Procedure 

An  iteration  procediare  was  used  to  obtain  the  ideal  impact  pres¬ 
sure,  p  ',  from  the  indicated  pressure,  pj.  An  initial  approximation 
was  made  equating  p^'  to  p^^.  Using  the  reservoir  stagnation  pressure, 
flow  properties,  such  as  Re 2,  r«  P2/P®.  M®.  P»>  and  were  obtained 

using  isentropic  flow  equations  along  with  the  value  for  p^'.  The  stag¬ 
nation  point  heat  transfer  was  obtained  using  Eqs.  (1-8),  (1-9),  and 
(I- 10).  The  Pyv  was  assumed  equal  to  p^,  thereby  permitting  the  solu- 
^on  K  from  Eq.  (1-7).  Equations  (1-3)  through  (1-6)  were  solved  for 
p  yielding  p^  which  was  substituted  in  Eq.  (1-7)  giving  a  new  k.  Equa¬ 
tions  (1-3)  through  (1-6)  were  repeatedly  solved  until  p^  did  not  change 
from  one  iteration  to  the  next.  Using  the  obtained  value  of  p^  in  Eqs. 
(I-l)  and  (1-2),  Pq  was  obtained  and  used  to  repeat  the  above  procedures 
until  the  value  of  Pq  was  constant  for  successive  iterations.  The  final 
value  of  Pq  is  the  ideal  impact  pressure  which  can  be  used  to  obtain  the 
corrected  free-stream  properties. 
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APPENDIX  II 

FLOW  DIAGNOSTICS  USING  EXCITATION  BY  AN  ELECTRON  BEAM 


THEORY  OF  THE  ELECTRON  BEAM  FLUORESCENCE  METHOD  OF  FLOW  DIAGNOSTICS 

In  -well-dociunented  papers  (Refs.  II- 1  and  II-2),  Muntz  has  pre¬ 
sented  in  detail  the  theory  underlying  the  excitation  of  a  low  density 
gas  by  an  electron  beam  and  use  of  the  induced  fluorescence  to  deter¬ 
mine  rotational  and  vibrational  temperatures  and  atomic  and  molecular 
specie  concentrations.  Therefore,  only  a  brief  summary  of  the  tech¬ 
nique  as  applied  to  helium  (He)  and  nitrogen  (N2)  is  presented  here. 

It  is  wen  known  that  the  interaction  of  an  electron  beam  of  energy 
of  approximately  20  kv  with  a  low  density  gas  can  be  characterized  by 
descriptions  of  two  general  coHisional  processes:  elastic  scattering 
and  inelastic  scattering.  Elastic  scattering  is  descriptive  of  the  devi¬ 
ation  of  the  path  of  the  beam  of  electrons  attributable  to  dose  encoun¬ 
ters  with  the  target  nudei  and  involves  no  excitation  of  the  gas  atoms. 

The  major  effect  of  dastic  scattering  is  the  spreading  of  the  beam, 
leading  to  the  necessity  of  the  determination  of  the  diameter  of  the  beam 
at  the  point  of  observation,  thereby  influencing  both  the  method  of 
current  measurement  as  well  as  the  technique  of  optical  detection.  The 
results  of  Center  (Ref.  11-3)  were  influential  in  selection  of  the  current 
measurement  method,  and  the  results  of  Maguire  (Ref.  II-4)  were  help¬ 
ful  in  selecting  parameters  for  the  optical  measurement  technique  for 
the  present  study.  Details  of  these  methods  of  measurement  are  dis¬ 
cussed  in  Section  3.  5  of  this  report. 

The  occurrence  of  inelastic  scattering  is  the  more  significant  pro¬ 
cess  of  the  two  for  the  purpose  of  flow  diagnostic  measurements,  for  as 
a  result  of  this  process  the  emission  of  radiation  results  which  can,  in 
principle,  be  related  to  gas  density  and  temperature.  Figure  II-l  shows 
a  sketch  of  two  possible  types  of  inelastic  collisions,  namely  ionizing 
and  nonionizing  collisions,  each  occurring  to  the  extent  determined  by 
its  partial  cross  section.  One  of  the  results  of  the  inelastic  collision 
can  be  excitation  of  the  atom  or  molecule  to  a  particular  level  of  an 
excited  state  from  which  it  decays  by  emission  of  radiant  energy^ 
intuitively,  one  expects  the  intensity  of  this  emission  to  be  related 
directly  to  the  gas  density,  thereby  affording  a  nonperturbative  optical 
method  of  density  measurement.  To  allow  use  of  such  a  method,  two 
major  criteria  must  be  met:  (1)  the  excitation  process  must  be  direct, 
i.  e. ,  free  from  secondary  production  sources  such  as  low  energy  elec¬ 
trons  and  cascading  transitions,  and  (2)  the  radiative  lifetime,  t^.. 
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Fig.  11-1  Schematics  of  Ionizing  and  Nonionizing  Collisions 
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must  be  sufficiently  short  that  effects  attributable  to  collisions  and 
flow  speed  are  ne^igible.  For  the  typical  hypersonic  flow  field 
environment  this  requires  to  be  less  than  approximately  lO”*^  sec, 


3^P 


Particularizing  the  discussion  to  He^  the  Hel  transition  of 
-►2^S  at  5016  A,  sketched  in  Fig.  n-2,  satisfies  both  these 
criteria  and  was  selected  for  use  in  this  test  for  the  purpose  of  per¬ 
forming  a  measurement  of  He  density.  It  is  now  necessary  to  cast 
into  quantitative  form  the  intuitive  feeling  of  the  relationship  of 
radiant  intensity  and  He  density.  For  sufficiently  low  gas  density,  for 
which  one  can  neglect  coUisional  effects,  it  is  shown  in  Ref.  11-2  that 
the  steady-state  intensity  per  unit  volume  is  given  by 


where 


and 


1(5016) 


Co  A{S016)  h  o  j  N(He) 
A<5016)  r,(3  *P) 


(n-1) 


A(5016)  =  Einstein  spontaneous  emission  coefficient 
h  =  Planck’s  constant 
c  =  Speed  of  light 

j  =  Electron  beam  current  density 

X(5016)  =  Wavelength  of  5016-Aline 

Tj,{3^P)  =  Radiative  lifetime  of  3^P  state 

Cq  =  Constant  independent  of  temperature  and 
density 

N{He)  =  Number  density  of  ground  state  ^Sq  He 
atoms. 


Equation  H-l  shows  the  intensity  I  to  be  directly  proportional  to 
the  number  density  of  the  ground-state  helium  atoms,  and,  therefore, 
a  measurement  of  I  will  yield  N(He)  if  by  calibration  the  multiplication 
factor  of  N(He)  is  determined.  If  this  is  done  Eq.  (H-l)  becomes 

1(5016)  =  CijN(He)  (II-2) 


where  is  ideally  independent  of  both  temperature  and  density. 
Details  of  the  calibration  required  to  determine  Cj  are  given  in  a 
following  section. 
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Referring  to  Fig.  II- 2,  it  is  shown  that  excited,  neutral  He  atoms 
in  the  level  not  only  radiatively  decay  by  the  transition 

3^Pl  -*•  2lSo  but  also  by  the  emission  of  an  ultraviolet  photon  of  wave¬ 
length  equal  to  53  7  A  in  the  transition  returning  the  atom 

to  the  electronic  groxmd  state.  Since  the  ultraviolet  (uv)  photon 
emitted  is  in  resonance  with  the  reverse  transition  S^Pj, 

there  exists  a  significant  absorption  coefficient  at  this  wavelength, 
once  again  returning  the  atom  to  the  S^P^  state  from  which  it  can  once 
more  radiatively  decay  by  emitting  a  photon  of  wavelength  of  either 
5016A  or  537 A.  This  has  the  effect  of  increasing  the  radiative  life¬ 
time  Tj,(3^P)  as  well  as  broadening  the  beam  by  diffusive  spreading. 
From  the  results  of  both  Maguire  (Ref.  n-4)  and  Phelps  (Ref.  II- 5), 
for  the  majority  of  density  values  encountered  in  the  test  the  effect 
was  to  broaden  the  beam  of  He  emission  species  to  approximately 
eight  "current  beam  diameters,"  and  allowance  was  made  for  this 
effect  by  adjustment  of  the  spectrometer  slit  height  to  ensure  complete 
imaging  of  the  beam. 


For  gas  density  values  such  that  the  collision  time  is  not  miich 
larger  than  t  one  must  alter  Eqs.  (II-l)  and  (II-2)  to  allow  for  collision 
decay  channels  for  excited  Hel  atoms;  this  is  shown  diagramraati- 
cally  in  Fig.  n-3  where  the  density-dependent  rate  constant  k  is  ex¬ 
pressed  in  units  of  cm^/sec.  For  the  simple  phenomenological  model 
of  that  figure  one  can  show  that  the  intensity-density  relationship  now 
becomes 


1(5016)  - 


r.2  j  N(iip) 

1  +  Ic  r,ulp)  N(llc) 


(n-3) 


where  C2  is  a  factor  independent  of  temperature  and,  because  of  the 
imprisonment  of  radiation,  is  dependent  on  the  number  density  N(He) 
to  a  slight  degree.  For  the  range  of  values  of  N{He)  encountered  dur¬ 
ing  this  test,  this  dependence  was  neglected.  It  has  been  assumed,  as 
before,  that  only  helium  is  present,  aind  the  more  general  case  of  a 
binary  mixture  is  considered  in  a  following  section. 


As  shown  by  the  results  presented  in  Fig.  11-13,  the  nonlinearity 
predicted  by  Eq.  (n-3)  becomes  significant  only  at  the  higher  density 
values  encountered.  For  values  of  N(He)  in  the  region  of  lO^^/cm^, 
Tj.(3^P)  is  known  to  be  approximately  7.  3  x  10“^  sec;  from  previous 
measurements  of  k  performed  locally  under  laboratory  conditions, 
the  quenchmg  constant  at  SOO'K  has  been  found  to  be  2.  0  x  10“^® 
cm^/sec.  From  Eq.  (n-3)  it  can  be  seen  that  at  300°K  quenching  be¬ 
comes  significant  for  values  of  N(He)  such  that  k*Tj.(3^Pj)  x  N(He)  =0.1 
or  for  N(He)  values  on  the  order  of  10^®/cm^.  If  one  assumes  hard 
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Fig.  11*2  Partial  Energy  Level  Diagram  for  Helium 


Deexcitation 


Fig.  11^  Phenomenoiogcal  Quenching  Radiative  Decay  Model 
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sphere  collision  theory  to  be  appropriate  for  the  description  of  the 
quenching  process,  it  is  found  that  k  varies  as  yjT,  and  is,  therefore, 
of  less  consequence  as  the  temperature  is  lowered  from  300°K  to 
those  values  encountered  in  the  test.  However,  it  is  difficult  at  this 
time  to  believe  that  k  behaves  in  this  manner  because  the  amount  of 
energy  transferred  by  the  quenching  process  -►  ^Sq  is  quite 
sizable,  as  Fig.  11-3  shows,  since  the  energy  separation  is  equivalent 
to  a  photon  with  a  wavelength  of  53  7A  which  is  approximately  23  ev; 
collisions  involving  such  a  large  transfer  of  energy  are  quite  inefficient. 
More  likely,  it  is  expected  that  a  quasi-resonance  transfer  is  the  mech¬ 
anism  responsible.  The  temperature  dependence  of  such  a  process 
varies  as  1/T,  becoming  more  important  as  the  temperature  is  lowered, 
but  only  further  experiments  will  decide  the  form  and  magnitude  of  the 
correction  to  be  applied,  because  such  data  are  not  in  hand  at  present. 
However,  one  may  appreciate  the  magnitude  of  the  possible  correction 
attributable  to  variation  of  k  with  temperature  if  one  assumes  k  to  vary 
as  1/T.  From  Eq.  (II-3)  it  can  be  seen  that  the  denominator  of  the 
right-hand  side  at  300“K  is  approximately  1.  14  at  N(He)  equal  to 

at  60TC,  for  the  same  number  density,  the  denominator  is 
approximately  equal  to  1.  70,  thereby  giving  approximately  a  30 -percent 
difference  in  density  values. 


EXCITATION  OF  GASEOUS  MOLECULAR  NITROGEN  BY  ELECTRONS 


Excitation  of  the  first  negative  system  of  N2'^  by  impact  of  high 
energy  electrons  upon  neutral  nitrogen  molecules  has  received 
extensive  treatment  in  the  literature  following  the  initial  report  of 
Muntz  (Ref.  II-l)  in  1961,  and  two  of  the  authors  of  this  report  have 
published  residts  of  an  investigation  of  the  appropriateness  of  the 
Muntz  model  of  excitation  as  well  as  an  extension  of  those  equations 
for  vibrational  temperature  measurements  at  higher  density  (Refs.  H-S 
and  7),  Reference  will  be  frequently  made  to  these  works  for  detailed 
derivations  of  the  equations  which  are  to  follow. 


Figure  n-4  shows  some  of  the  energy  levels  of  N2  and  N2^  mole¬ 
cules  along  with  the  designation  of  the  radiative  molecular  band  sys¬ 
tem  denoted  as  the  N2'*'  first  negative,  or  N2’*'  (l-)y'y'o  in  which  the 

subscripts  v'v"  are  used  to  designate  the  upper  and  lower  vibrational 
states,  respectively.  Shown  in  Fig.  H-S  is  an  enlarged  diagram  of  the 
excitation  emission  process  caused  by  the  inelastic  collision  between 
a  primary  electron  ep  and  a  neutral  nitrogen  moleculej  i.  e. 
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N2  X^2(v 


N2^B22(v',K') 
- -  N2'^X*S(v 


+  ep 

",K") 


+  Cs  -  AE 
+  hv 


(n-4) 


14-1 


N2A3e 


AEDC-TR-71-25 


N2+b2c 


N2X  £ 


K"'=  0;  V- 


Fig.  H-5  Excitation-Emission  Diagram 


143 


AEDC-TR-71-25 


where  K  denotes  the  rotational  quantum  number  of  the  molecule,  and  it 
is  to  be  noted  that  dipole  excitations,  =  ±  1,  are  assumed  to  be  the 
production  mechanism  and  dipole  emission  the  process  of  radiative 
decay.  With  these  assumptions,  and  further  assuming  the  negligibility 
of  collisional  deexcitation  and  excitation  by  the  secondary  electrons, 
designated  in  Eq,  (II-4)  by  eg,  one  finds  the  steady-state  intensity 

v'K' 

I  .  of  an  emission  line  for  the  R  branch  transition 
N2+b=2:(v;k")  N2+x*2(v"K'  =  k"-  i) 


'  r  V  V  ^ 

'■k  -I  -1  CjN{N2)  "(K'+Dhc/kTH 

'[G],'=o  ■  (p.r.)^  ■  {n-5) 


where 

C  =  constant 

N(N2)  =  total  N2X^21  number  density 

Bq  =  rotational  constant  for  the  v'"  =  0  state  of  N2 
k  “  Boltzmann's  constant 
Tj^  =  rotational  temperature 

(p.f.  )/j«  =  total  partition  function  of  N2X^E  which  possesses  a  vibra¬ 
tional  temperature  and  a  rotational  temperature  Tp 

v'K' 

A  ,  =  wavelength  of  the  emission  line 

V  -i 


[Glo 


(K'+  +  k'  e+2^«Kkc/kTR 

2K'+  i 


and  [G]q  is  known  as  the  Mimtz  G  factor,  or  simply  as  the  G  factor. 
Taking  the  natural  logarithm  of  both  sides  one  finds 


which  shows  that  the  experimental  intensity  data,  when  cast  in  the  form 
of  the  left-hand  side  of  the  equation  and  plotted  versus  K'(K*  +  1)  in 
principle  yield  a  straight  line  whose  slope  equals  -Bohc/h  Tp,  thereby 
giving  Tj^.  Since  [g]q  contains  Tj^,  Eq.  (II-3)  must  be  solved  by  itera¬ 
tion,  but  in  practice  convergence  is  rapid. 
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Performing  a  summation,  either  mathematically  or,  in  practice, 
electronically,  of  the  rotational  line  intensity  values  over  the  vibra¬ 
tional  band  v'v”  one  obtains  the  vibrational  band  intensity 


I  =  1  ^  = 

T  TV 


Cl  N(N2)  j  p(v'v")  2 

V  =Q _ 

(p.f.)v  '"X*  %  -  £  \  "1 

T  ''■^=0 


(II- 7) 


where 

(p,  f.  )y#f#  =  vibrational  position  function  of  the  N2 
molecule 

”  average  wavelength  of  the  vibrational  band 
G(v'")  =  vibrational  constant  of  the  N2X^E  molecule 
p(v'v")  =  vibrational  band  strength  for  emission 
p(vV"')  =  vibrational  band  strength  for  excitation 


If  one  now  obtains  the  ratio  of  two  vibrational  bands  originating 
from  different  upper-state  vibrational  levels  i  and  k,  respectively, 
and  terminating  in  levels  j  and  1,  respectively,  one  finds 


!il 

Ikl 


p(t  1) 


i-n 


2  p(ka)  r,(k) 

•=B 


(n-8) 


where  r^Ci)  =  1/  E  p(i  s)/ g  =  radiative  lifetime  of  level  i. 
s=0 


It  is  well  known  that  the  vibrational  band  strength  p  (i  j),  for 
example,  can  be  written  as 

p(ij)  =  Re^  (rij)  q(ij)  {n-9) 

where  q(i  j)  is  the  calculable  Frandc- Condon  or  overlap  factor  for  the 
i  j  transition  and  is  given  in  terms  of  the  vibrational  configurations 

by 

q(ij)  =  l<^ij>l"  (n-10) 

where  |i  >  denotes  the  appropriate  vibrational  eigenket  vector. 
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AlSOj 

Re2(fip  =  |<i|  R^(r)  |j>!2  (II- 11) 

Z 

where  electronic  transition  moment,  The  terms 

|t)>  and  (d>  represent  the  final  and  initial  electronic  state  eigenkets. 
respectivelyj  and  r^g  is  the  displacement  vector  relative  to  the  origin  of 
the  ith  electron. 

The  electronic  transition  moment  is  calculable  only  for  the  sim¬ 
plest  molecules  and,  as  a  restilt,  must  be  experimentally  determined 
for  the  N2^(1-)  system. 

The  r- centroid  r^j  is  defined  by 

fij  =  <ilrlj>,<ilj>  (11-12) 

and  represents  an  ’’average"  internuclear  distance  r  for  the  transition. 

Table  II- 1  lists  the  calculated  values  of  q(ij)  and  r^j  as  well  as 
experimentally  determined  values  of  p(ij).  It  is  noted  that  for  three 
bands  for  which  i-j  =  k-1  «=  Av  and  i  k,  the  values  of  the  first 
r -centroids  are  very  nearly  equal,  as  are  the  wavelengths  Xij  and 
Xjjj.  Since  the  r-centroids  are  very  nearly  equal,  then  RgCfy)  a 

Re(f‘ki)  and 

p(ij>  _  qUj) 
p{kl)  q(kl) 

showing  that  use  of  the  two  bands  ij  and  kl  has  the  results  of  minimiz¬ 
ing  the  inaccuracy  of  a  vibrational  temperature  measxirement,  since 
Rg(f^j)  usually  has  large  errors  associated  with  it.  Additionally,  since 
X-ij  S£  Xjj]^  for  i-j  =  k  -  1 ,  selection  of  these  two  band  systems  will  mini¬ 
mize  the  error  associated  with  making  spectral  sensitivity  corrections. 
For  these  two  reasons,  the  band  pairs  (01),  (12),  and  (02),  (13)  were 
used  for  vibrational  temperature  measurements  using  the  N2'*’  (1-) 
system. 

One  should  note  that  Eq.  (II-8)  contains  not  only  emission  band 
strengths,  pjj,  but  also  excitation  band  strengths,  p  (is),  which  have 
the  same  mathematical  basis  as  the  emission  band  strengths;  unfor¬ 
tunately,  no  work  has  been  performed  to  determine  their  values.  Be¬ 
cause  of  this  ignorance,  and  with  no  other  justification,  excitation 
Franck-Condon  factors,  q(is),  are  used  in  their  stead,  and  these 
parameters  are  listed  in  Table  11-2. 
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v' 

v" 

q{v'v"  ) 

p(v'.v"  ) 

r(v'v" 

0 

0 

0. 6570 

0. 6364 

1.0995 

0 

1 

0. 2982 

0. 2611 

1.  1523 

1 

0 

0. 2546 

0.  2693 

1. 0512 

1 

1 

0.2316 

0.  2189 

1.  1124 

2 

0 

0.  0686 

0.  0776 

1. 0227 

2 

1 

0.  2838 

0.  2971 

1. 0572 

3 

0 

0. 0157 

0.  0189 

0.  9543 

3 

1 

0. 1299 

0.  1458 

1.  0177 

4 

0 

0. 0033 

0. 0004 

0.  8923 

4 

1 

0.0418 

0.  0498 

0. 9832 

TABLE  11-2 

EXCITATION  FRANCK-CONDON  FACTORS  q(v',v"  ) 


qW'v"  ) 


0 

0 

1 

1 

2 

2 

3 

3 


0 

1 

0 

1 

0 

1 

0 

1 


0. 8794 
0. 1184 
0.  1063 
0.  6874 
0.  0128 
0. 1606 
0.  0014 
0. 0290 
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The  radiative  lifetimes  7^(0)  and  Tp(  1)  are  required  to  evaluate 
Eq.  (II-8),  and  for  this  purpose  the  value  of  Tj.(0)  =  6.  58  x  10"^  re¬ 
ported  by  Bennett  and  Dalby  (Ref.  II-8)  is  used.  Using  results  of 
experimental  intensity  data  Tj.(1)  was  found  to  be  6.50  x  10"^  sec. 
With  these  values  the  ratio  of  intensities  of  the  two  just -mentioned 
band  pairs  are  shown  in  Fig.  II-6  as  a  function  of  Ty;  with  these 
curves  measured  ratios  of  band  intensity  can  be  related  directly  to 
the  vibrational  temperature. 


Rettirning  to  Eq.  (IE- 7),  one  sees  that  the  band  intensity  can 
be  written  as 


1 


C2  j  N(N2) 


'^0  * 
V  =0 


hc/kx. 


(n-13) 


where  C2  is  a  constant,  thereby  enabling  one  to  make  a  determination 
of  the  neutral  nitrogen  number  density  N(N2)  by  simply  measuring  the 
integrated  band  intensity,  so  long  as  Ty  is  known. 


Selecting  the  most  intense  vibrational  band  y'  =  v"  »  0  for  use  in 
Eq.  (n-13)  and  the  measurement  of  N(N2),  Fig.  II-7,  shows  Iqo  sis  a 
function  of  N(N2)  for  various  values  of  Ty,  thereby  demonstrating  not 
only  the  linear  dependence  of  Iqo  on  N(N2)  but  also  the  errors  to  be 
encountered  by  an  ignorance  of  Ty. 

In  the  event  the  density  N(N2)  is  sufficiently  large  that  kTpN(N2) 
becomes  significant  relative  to  unity,  as  with  He,  effects  attributable 
to  addition  of  collision  decay  channels  must  be  considered.  Figure  11-8 
shows  a  phenomenological  model  of  the  radiative  and  collisional  deexci¬ 
tation  processes  for  the  vibrational  levels  0  and  1  of  the  N2*^  E  state 
where  kQ(N2'*',  N2)  and  k^  (N2'*',  N2)  represent  the  density-dependent 
rate  constants  for  N2  collision  quenching  of  the  v  =  0  and  1  levels, 
respectively.  With  this  model  one  can  easily  show  that  the  intensity  of 
the  v’  =0,  v"  =0  band  can  be  written  as 

=  I„“  .  (N(N2)/[1  +  ko  (N2MV2)  r  (0)  N(N2)]  (n-14) 

where  the  vibrational -temperature-dependent  quantity  1°^  is  the  limit 
of  as  N(N2)  approaches  zero. 

Values  of  kQ(N2'*^,  N2)  and  ki(N2'*’,  N2)  have  been  estimated  to  be 
2. 1  X  10"^^  cm3/sec  and  3. 1  x  10-10  cm^/sec,  respectively,  at  300'K. 
With  the  same  qualifications  for  the  temperature  dependence  of  the 
quenching  constants  as  stated  for  He,  one  can  now  employ  Eq.  (11-14) 
for  a  measurement  of  N(N2)  for  those  regions  of  N(N2)  for  which 
collisional  effects  are  important,  the  upper  N(N2)  lim.it  being  deter¬ 
mined  by  the  decreasing  slope  of  Iqq  versus  N(N2)  and  the  slope  finally 
approaching  zero  at  several  torr  at  room  temperature. 
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To  conclude  remarks  concerning  pure  nitrogen,  the  effects  of 
collisions  upon  the  vibrational  temperature  determination  are  obtained 
by  generalizing  Eq.  (11-8)  using  Eq.  (n-14)  to  yield 

loi  _  /loj  \  ^  [l  +ki(N2^  N2>r,(L)N(N2)l 

fll  ~  Vll  /  [l  +k<,(N2^  N2>M0)N(N2)]  (n-15) 


where 


is  the  low  density  limit  of  Eq. 


(n-15)  and  is  identical  with 


the  right-hand  side  of  Eq,  (n-8).  It  is  emphasized  that  dependence  of 
Eq.  (n-15)  upon  the  kinetic  temperature,  T^,  is  at  present  unknown, 
and  when  employed,  it  must  be  done  with  caution. 


EXCITATION  OF  A  BINARY  MIXTURE  OF  HELIUM  AND  NITROGEN 


In  the  event  the  gas  sample  is  a  binary  mixture  of  He  and  N2,  the 
preceding  equations,  notably  Eqs.  (n-2),  (n-6),  (n-8),  and  (n-13),  are 
valid  as  long  as  coUisional  effects  are  negligible.  Over  the  density 
regime  for  which  coUisional  effects  first  become  discernible,  one 
expects  Eqs.  (U-S),  (11-14),  and  (11-15)  to  serve  adequately  as  first- 
order  approximations.  However,  for  sufficiently  high  values  of  density 
the  above  equations  must  be  replaced  by  expressions  in  accord  with  the 
approximate  model  as  shown  in  Figs,  n-9  and  H-IO.  These  appropriate 
equations  are 


1(5016)  = 


C2  •  j  •  N(He) 

1  +  [k(He,  H«)  N(He)  +  k(He,  N2)  N(N2)]  r(3>  P) 


(n-16) 


loo  =  loo"  N(N2)/I1  +  1o(N2^  N2)N(N2)  +  ko(N2^  He)N(He)l  ■  r^CO) 


'OJ 


hi 


il  ■^■[kl(N2^  Ng)  N(N2)  +  ki(N2^  He)  N(He)]  ^1) 
ll  +  [k„<N2^  Ng)  N(N2)  +ki(N2'^,  He)  N(IIe)l  r,(0) 


(n-i7) 

(n-18) 


Of  the  quenching  constants  found  in  the  preceding  three  equations. 

He),  ki(N2'*',  He),  and  k(He,  N2)  are  unknown  at  present  and 
can  only  be  estimated  if  needed.  Table  U-S  lists  the  values  of  those 
rate  constants  that  were  previously  mentioned  as  having  been  obtained 
at  SOOTK  and  also  shows  values  of  those  unknown  quenching  rates  that 
are  estimated  using  elastic  scattering  collision  rates  of  charged- 
uncharged  pairs  of  species  (Ref.  II- 9)- 


Rewriting  Eq.  (11-18)  as 


hi 


■fej 


R 


(U-19) 
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a.  Phenomenological  Model  of  Collision 


He  or  N2 


(He)=*' 


Before 


Lower  Level(s) 


He  or  No 

-  ^ 

(Hel)* - 

After 


b.  Schematic  of  Collision  Process 
Fig.  11-9  Collisional-Radiative  Deexcitation  Model  of  Hel 
(3’  Pi )  in  Mixture  of  Nitrogen  and  Helium 
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He  or  N2  N2'^  He  or  N2 


b.  Schematic  of  Collision  Process 
Fig.  Ii-10  Collisional-Radiative  Decay  Scheme  of  N2'^B2S 
in  a  Mixture  of  Nitrogen  and  Helium 
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TABLE  11-3 

QUENCHING  RATE  CONSTANTS 


Rate  Constant 

T.  “K 

Magnitude,  cm^/sec 

k(He,  He) 

300 

1.98  X  10“^® 

ko(N2+.  Ng) 

300 

2.  1  X  10"^° 

ki(N2'".  N2) 

300 

3. 1  X  10"^° 

k(He.  N2'^) 

(10-10)* 

He) 

(lO'lO)* 

ki(N2‘^,  He) 

(  10"^°)* 

*  Estimated 
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wnere  R  is  the  density -dependent  factor  of  Eq.  (11-18),  one  sees  that, 
using  the  estimated  rate  constants  of  Table  II- 3 

1  +  0.02015  NCNa)  +  0.0065  N(He) 

1  +  0.01382  NlNa)  +  0.0066  N{He)  (11-20) 

H  c  o 

where  the  number  densities  are  expressed  in  units  of  10  /cm'^. 

Figure  11- 11  shows  R  as  a  function  of  the  nitrogen  number  density 
N(N2),  assuming  N(He)  is  zero.  For  the  range  of  N(He)  values  en¬ 
countered  in  the  test,  the  correction  attributable  to  quenching  collisions 
with  helium  atoms  is  negligible — a  result  which  was  ensured  by  assum¬ 
ing  the  equality  of  values  of  kQ(N2''',  He)  and  ki(N2'*',  He). 

The  density  correction  for  the  rotational  temperature  data  has  not 
been  discussed  because  none  of  the  quenching  constants  for  the  individ¬ 
ual  rotational  levels  are  at  present  known.  If  sufficient  interest  and 
requirements  exist,  this  determination  is  straightforward. 

CALIBRATION  OF  ELECTRON  BEAM 
HELIUM  NUMBER  DENSITY  CALIBRATION 

Calibration  of  the  device  consists  of  determining  the  constant  C  in 
Eq.  (n-2)  which  relates  the  intensity  of  the  501 6-A  line  to  the  current 
density  and  helium  number  density.  The  calibration  was  performed 
statically  in  the  tunnel  using  the  system  shown  in  Fig.  n-12.  The  sys¬ 
tem  consists  of  an  electron  drift  tube  extension,  test  chamber,  vacuum 
pumps  and  valves,  bleed  valve,  helium  supply,  Equibar  pressure  meter, 
and  retractable  beam  receiver  cup.  The  drift  tube  and  test  chamber 
were  bolted  to  the  orifice  plate  of  the  electron  gun  system.  Test  cham¬ 
ber  pressure  was  measured  with  the  Equibar  pressure  meter  (a  portable, 
differential,  variable-capacitance  micromanometer).  Reference  pres¬ 
sure  for  the  pressure  meter  was  the  electron  gun  section  pressure 
(less  than  1  x  10"^  torr).  Beam  current  was  read  with  a  digital  volt¬ 
meter. 

The  calibration  system  (with  heliiun  supply  valve  closed)  was  first 
evacuated  to  a  pressure  of  approximately  0,  2  torr.  The  bleed  valve 
was  then  closed,  and  the  helium  supply  valve  opened.  A  pressure  of 
40  psia  was  then  set  on  the  high  pressure  side  of  the  bleed  valve.  The 
vent  valve  on  the  high  pressure  side  of  the  bleed  valve  was  opened  for 
approximately  1  min  to  further  purify  the  helium  to  be  leaked  into  the 
test  chamber.  With  the  bleed  valve  dosed,  the  remainder  of  the  cali¬ 
bration  system  was  pumped  down  to  less  than  5  x  10“^  torr  with  a  leak  . 
rate  of  less  than  1  x  10~^  torr/ min.  Helium  was  then  leaked  into  the 
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Fig.  11-12  Schematic  Drawing  of  Tunnel  M  Electron  Beam  Installation  and 
Helium  Calibration  S/stem 
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test  chamber  through  the  bleed  valve  until  the  desired  test  chamber 
pressure  was  achieved.  The  electron  beam  was  then  injected  into  the 
chamber.  The  intensily  of  the  5016-Aline  was  measured  using  the 
integration  technique  described  in  Section  3.  5.  For  each  pressure 
setting  the  line  intensity  was  measured  twice,  and  the  beam  current 
was  measured  with  the  retractable  receiver  cup  positioned  at  the  point 
along  the  electron  beam  that  was  focussed  on  the  spectrometer  slit. 
There  was  a  small  difference  in  beam  area  distribution  between  the 
test  chamber  and  actual  plume  conditions.  However,  from  the  data  of 
Center  (Ref.  11-3),  for  a  number  density  of  1.  5  x  lO^^/cm^,  the  beam 
diameter  is  approximately  15  mm.  Since  the  magnification  factor  for 
the  optical  system  was  1/2,  this  gave  an  effective  diameter  of  7.  5  mm 
which  was  easily  collected  by  the  20-mm-high  spectrometer  slit. 

Two  helium  calibrations  were  made  and  are  shown  in  Fig.  11-13. 
The  calibrations  are  plotted  as  the  ratio  of  5016-A  intensity  to  beam 
current  versus  helium  pressure  at  300'IC.  This  ratio  is  a  fairly  linear 
function  of  helium  density  as  it  should  be  over  the  calibration  density 
range  in  Fig.  11-13. 

NITROGEN  NUMBER  DENSITY  CALIBRATION 

Nitrogen  number  densities  are  related  to  the  intensity  through 
Eq.  (n-13).  The  reference  density  was  taken  as  the  Tunnel  M  center- 
line  free-stream  density.  The  N2'*’[l-]  (0,0)  band  intensity  was  meas¬ 
ured  using  the  band  integration  technique,  and  beam  current  was  meas¬ 
ured  at  both  the  auxiliary  and  flat-plate  positions.  With  each  measure¬ 
ment  of  band  intensity  and  current,  readings  of  the  ttmnel  conditions 
were  made.  Normally  three  of  four  band  measurements  were  made  to 
give  an  average  ratio  of  band  intensity  to  beam  current,  and  the  tunnel 
data  reduction  program  provided  a  calculation  of  nitrogen  nimiber 
density.  With  the  valid  assumption  that  the  relation  of  band  intensity 
to  N2  number  density  was  linear  for  the  number  density  range  in  this 
test,  this  single-point  calibration  technique  provided  the  necessary 
calibration  of  band  intensity  versus  number  density.  This  calibration 
procedure  was  repeated  each  time  that  nitrogen  number  density  meas¬ 
urements  were  made  in  the  plirnie.  Slit  settings,  photomultiplier 
supply  voltage,  and  lock-in  amplifier  scales  were  kept  the  same  during 
tunnel  measurements  as  they  were  during  the  calibration. 

SPECTRAL  SENSITIVITY  CALIBRATION 

In  order  to  check  the  relative  spectral  sensitivity  of  the  optical 
system,  a  calibrated,  tungsten-ribbon,  standard  lamp  was'  set  inside 
Tunnel  M.  With  the  lamp  turned  on,  the  test  section  was  closed  and 
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Relative  Senaltivlty 


Fig.  11-14  Relative  Sensitivity  of  Tunnel  M  Optical  System 
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the  lamp  ribbon  was  focussed  onto  the  slit  of  the  1-meter  spectrometer. 
In  this  manner  the  relative  spectral  sensitivity  of  the  entire  optical  sys¬ 
tem  was  checked.  Two  wavelength  scans  from  3500  to  6500  A  were 
made,  and  the  relative  sensitivity  of  the  Tunnel  M  optical  system  is 
shown  in  Fig.  11-14. 

As  explained  in  Section  4.4  a  ratio  of  band  intensities  was  used  for 
a  determination  of  vibrational  temperature.  As  can  be  seen  from 
Pig.  11-14  a  small  relative  sensitivity  correction  had  to  be  applied  to 
these  ratios.  For  the  rotational  scans  of  the  R  branch  of  the 
^2**^  &"]  (0,0)  band  only  about  a  1-percent  relative  sensitivity  varia¬ 
tion  is  evident,  and  therefore  neglected. 
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